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ABSTRACT 


Linear  equations  of  motion  with  the  hydrodynamic 
coefficients  approximated  by  strip  method  are  used  to  obtain 
the  motion  of  small-waterplane-area,  twin-hull  (SWATH)  ships 
in  regular  and  irregular  waves.  The  couplings  among  different 
modes  of  motion  are  separated  into  three  independent  groups, 
surge,  heave-pitch  and  sway-roll-yaw,  ar.d  the  resulting  equations 
of  motions  are  treated  in  the  frequency  domain. 

The  hydrodynamic  coefficients  of  predominantly  nonviscous 
nature,  added  mass  and  wavemaking  damping,  are  obtained  by 
strip  theory.  The  effects  of  wave  diffraction,  viscous  damping, 
and  stationary  stabilizing  fins  are  included  in  the  equations  of 
motion  within  the  frame  work  of  the  linear  frequency  response 
of  a body  to  waves.  The  assumption  involved  in  the  evaluation 
of  various  hydrodynamic  coefficients  in  the  equation  of  motion 
are  specified. 

The  validity  of  the  theoretical  prediction  of  motion  is 
checked  by  correlating  with  the  existing  model  experimental 
results,  and  the  relevant  discussions  are  made.  In  general,  the 
theoretical  prediction  is  found  satisfactory  except  the  case  of 
near  zero-encounter  frequency  which  can  occur  when  a ship 
cruises  with  certain  speed  in  stern-quartering  waves.  A future 
effort  for  improvement  of  the  theoretical  prediction  of  SWATH 
ship  motion  in  stern-quartering  waves  is  recommended. 


ADMINISTRATIVE  INFORMATION 

This  study  was  sponsored  by  the  Naval  Sea  Systems  Command  as  part  of  the  High- 
Performance  Vehicle  Hydrodynamic  Program  of  the  Ship  Performance  Department,  David  W. 
Taylor  Naval  Ship  Research  and  Development  Center  (DTNSRDC).  Funding  was  provided 
under  Task  Area  SF  4342 1 202,  Task  18247,  Work  Unit  1507-200. 

INTRODUCTION 

A small-waterplane-area,  twin-hull  (SWATH)  ship  consists  of  submerged  twin  hulls,  an 
above-water  hull,  and  connecting  struts.  Figure  1 is  schematic  of  a typical  SWATH 
configuration. 
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TOP  VIEW 


Figure  1 - Configuration  of  a SWATH  Ship 
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Ihc  concept  ot  the  SWATH  evolved  through  many  years  and  its  development  is  well 
documented.'  5 One  ol  the  major  advantages  of  the  SWAIH  concept  lies  in  its  improved 
seakeeping  quality  compared  to  a monohull  ship  of  equivalent  displacement.  The  improve- 
ment is  contributed  mainly  by  the  larger  natural  periods  of  heave  and  pitch  modes  due  to  the 
small  waterpiane  area  and  by  the  smaller  wave-exciting  forces  and  moments  due  to  the  semi- 
submersible  configuration.4  The  larger  natural  periods  reduce  the  likelihood  that  resonant 
motions  will  be  excited  in  moderate  seas,  and  the  smaller  wave-exciting  forces  and  moments 
not  only  reduce  the  motion  but  also  provide  a greater  possibility  for  controlling  its  effects 
through  the  use  of  fins. 

I he  overall  assessment  of  a new  hull  concept  versus  existing  ships  is  highly  complicated 
and  definitely  beyond  the  scope  of  the  present  work.  The  aim  here  is  to  provide  an  analytical 
tool  which  can  be  utilized  to  assess  certain  aspects  of  the  seakeeping  qualities  of  SWATH  ships 
and  which  may  later  be  used  to  design  specific  hull  forms. 

Several  analytical  methods  have  been  developed  in  various  parts  of  the  world  in  the  past 
few  years  for  predicting  the  motions  of  monohulls.  Most  of  them  use  the  so-called  strip 
theory  to  evaluate  the  hydrodynamic  coefficients  associated  with  motion  of  ships.  The  strip 
theory  is  based  on  the  assumption  of  two-dimensionality  of  fluid  motion  surrounding  a cross 
section  ol  a ship.  A typical  application  of  the  strip  theory  is  well  described  in  Salvesen  et  al.5 

Equations  of  motion  and  strip  theory  similar  to  those  used  in  the  cited  reference  5 are 
employed  here.  However,  the  present  work  requires  a considerable  number  of  additional 
hydrodynamic  coefficients  which  are  normally  absent  in  the  equations  of  motion  for  mono- 
hull ships.  Three  distinct  considerations  should  be  included  in  the  hydrodynamic  coefficients 
for  SWATH  ships: 

1.  Hydrodynamic  interactions  between  the  two  hulls. 

2.  The  viscous  damping  effects.  This  cannot  be  neglected  because  for  SWATH  ships  it  is  of 
the  same  order  of  magnitude  as  wavemaking  damping.  As  expected  for  a lightly  damped 


1 Lang,  T.G.  and  D.T.  Higdon,  “S3  Semi-Submerged  Ship  Concept  and  Dynamic  Characteristics,'’  AIAA/ 
SNAME/USN  Advanced  Marine  Vehicles  Meeting,  Annapolis.  Maryland  (Jul  1072). 

A complete  listing  of  references  is  given  on  pages  78-80. 

'Leopold.  R.  et  al„  “The  Low  Water  Plane  Multi-Hull  Principles.  Status,  and  Plans."  AIAA/SNAME/USN 
Advanced  Marine  Vehicles  Meeting,  Annapolis,  Maryland  (Jul  1072). 

^Hawkins,  S.  and  T.  Sarchin,  “The  Small  Waterplane-Area  Twin  Hull  (SWATH)  Program  - A Status  Report,” 
AIAA/SNAME  Advanced  Marine  Vehicles  Meeting,  San  Diego,  California  (Jul  1074). 

4Motora,  S.  and  T.  Koyama,  “Wave-Excitationless  Ship  Forms,"  6th  Naval  Hydrodynamic  Symposium. 
Washington,  D.C.;  proceedings  published  by  the  Office  of  Naval  Research,  pp.  383-41  I (1066). 

5Salvesen.  N.  et  a I.,  “Ship  Motion  and  Sea  Loads.”  Trans.  SNAME,  Vol.  78  pp.  250-287  (1070). 
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system,  neglect  of  viscous  damping  effects  would  yield  unrealistically  large  motion  amplitudes 
at  the  resonant  frequencies. 

3.  Effect  of  stabilizing  fins.  A SWATH  can  become  unstable  in  the  vertical  plane  at  and 
beyond  a certain  speed  because  of  the  small  waterplane  area.  This  is  mainly  caused  by  a 
destabilizing  pitch  moment,  often  referred  to  as  the  Munk  moment,  which  is  proportional  to 
the  square  of  the  forward  velocity  of  the  body.  Hence,  stabilizing  fins  are  necessary  to 
augment  the  stability  of  the  ship.  The  stabilizing  fins  can  also  provide  much  needed  damping 
effects  to  the  motion  of  a SWATH  through  the  hydrodynamic  lift  generated  by  an  angle  of 
attack  which  results  from  the  combination  of  the  forward  motion  and  the  vertical  motion  of 
the  tins.  In  order  to  satisfactorily  compute  the  motions,  the  hydrodynamic  effects  contributed 
by  the  fins  should  be  included  in  the  equations  of  motion. 

A theoretical  prediction  of  motion  and  sea  loads  of  twin-hull  ships  has  been  presented  by 
Pien  and  Lee.6  Since  that  time,  significant  improvements  has  been  achieved  in  the  prediction 
of  SWATH  ship  motions.  These  mainly  involve  (I ) a more  adequate  representation  of  the 
viscous  damping  in  the  equations  of  motion  and  (2)  capability  to  predict  the  motion  of  a 
SWATH  with  stationary,  horizontal  fins. 

The  present  report  provides  a detailed  description  of  the  theoretical  analysis  of  the 
motion  of  SWATH  ships  in  waves  and  the  incorporation  of  that  description  into  computer 
programs.7  The  first  section  of  the  report  describes  the  equations  of  motion,  their  solutions, 
and  the  statistical  averages  of  motion  in  irregular  waves.  The  second  section  presents  the 
derivation  of  the  hydrodynamic  coefficients.  The  third  presents  a correlation  of  theoretical 
results  with  model  experimental  results  and  pertinent  discussions  of  these  results. 

Although  comparisons  of  computed  and  measured  values  of  the  hydrodynamic  coefficients 
and  the  motion  suggest  that  the  present  moiion  prediction  method  is  in  general  satisfactory, 
it  is  felt  that  further  investigations  should  be  pursued  on  such  aspects  as  viscous  damping 
effects,  fin-body  hydrodynamic  interactions  under  a free  surface,  and  prediction  of  motion  in 
stern-quartering  seas. 


6Picn.  P.C.  and  C.M.  Lee.  “Motion  and  Resistance  of  a Low-Waterplane  Catamaran,”  9th  Naval  Hydrodynamic 
Symposium,  Paris,  France,  proceedings  published  by  the  Office  of  Naval  Research,  pp.  463-545  (1972). 

7McCreight,  K.K.  and  C.M.  Lee,  “Manual  for  Mono-Hull  or  Twin-Hull  Ship  Motion  Prediction  Computer 
Program,”  DTNSRDC  Report  SPD-686-02  (1976). 
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EQUATIONS  OF  MOTION  AND  SOLUTIONS 


The  formulation  of  the  equations  of  motion  will  be  limited  to  linear  rigid  body  dynamic 
responses  of  the  body  to  harmonic  exciting  forces  or  moments.  The  exciting  forces  and 
moments  are  assumed  to  be  solely  contributed  by  free-surface  waves;  to  justify  the  linear 
response  of  the  body,  the  wave  slopes  are  assumed  to  be  small.  The  ocean  in  which  the  ship 
is  underway  is  infinitely  deep  and  has  no  appreciable  currents  or  winds  which  would  invalidate 
the  linear  response  assumption.  With  such  an  ocean  environment,  it  can  be  safely  assumed 
that  the  ship  can  maintain  a constant  mean  speed  and  a straight  mean  course.  The  submerged 
portion  of  the  ship  hull  is  assumed  to  be  sufficiently  slender  that  the  rate  of  lengthwise 
variation  of  the  hydrodynamic  pressures  is  small  compared  to  the  rate  of  the  variation  in  the 
girthwise  direction  in  a cross  section  of  the  ship. 

The  reference  frame  for  which  the  equations  of  motion  are  to  be  formulated  is  a right- 
handed  Cartesian  coordinate  system  Oxyz  which  translates  on  the  mean  path  of  the  ship  with 
the  ship  speed.  The  origin  is  located  on  the  undisturbed  free  surface,  and  the  Oxy  plane 
coincides  with  this  surface  When  the  ship  is  at  its  mean  position,  the  Oxz-plane  contains  the 
longitudinal  plane  of  symmetry  of  the  ship,  the  positive  Ox-axis  is  directed  toward  the  bow. 
and  the  positive  Oz-axis  is  directed  vertically  upward.  Figure  l shows  the  coordinate  system. 

In  general,  the  linearly  coupled  motion  of  a ship  can  be  expressed  in  the  form 

6 

Z t(Mik  + Aik » *k  + Bik  h + Cik  *k  1 = Fi<e)  ' ' > 

k = 1 

Here  i is  1 for  surge.  2 for  sway.  3 for  heave.  4 for  roll.  5 for  pitch,  and  6 for  yaw;  £k  is  the 
linear  or  angular  displacement  of  the  ship  from  its  mean  position  in  the  kth  mode.  |k  and  £k 
are  respectively  the  acceleration  and  velocity.  Mjk  is  the  mass  matrix.  Ajk  is  the  added  inertia. 
Bjk  is  the  damping.  Cjk  is  the  restoring  or  spring  constant.  F^e)  is  the  complex  amplitude  of 
the  wave-exciting  force,  j is  the  imaginary  unit  associated  only  with  a harmonic-time  function, 
and  oj  is  the  wave-encounteT  frequency.  The  expression  added  mass  (or  inertia)  which  will  be 
frequently  employed  in  this  report  refers  to  the  hydrodynamic  coefficient  associated  with  the 
acceleration  of  the  body  so  that  the  product  of  the  two  provides  a hydrodynamic  equivalent 
of  the  inertia  force  or  moment  of  the  body.  To  be  compatible  with  the  complex  expression 
on  the  right-hand  side  of  Equation  (I ).  the  motion  displacements  £k  are  assumed  to  be 
complex  functions  given  by 

$k(t>  = $koe'jWts<$kc+j  We'jW‘  <-> 
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where  £kc  and  £ks  are  real  functions  and  it  is  understood  that  whenever  a product  involved 
with  e'f001  appears,  only  the  real  part  of  the  product  would  be  recognized,  i.e., 

£k(t)  = £kc  cos  cJt  + £ks  sin  cot 

The  surge  mode  is  assumed  to  be  decoupled  from  the  rest  of  the  modes,  and,  furthermore, 
the  added  mass,  damping,  and  diffracted  wave  force  in  the  x-direction  arc  assumed  to  be  negligi- 
bly small.  Thus,  the  equation  of  surge  motion  can  be  expressed  by 

M = Fj(FK)  e'jWI  (3) 

where 

M = Mjj,  i = 1,  2,  3 

is  the  mass  of  the  body  and  appendages,  if  any.  Substitution  of  Equation  (2)  into  (3)  yields 

p (FK) 

^1°  ~ 2 M (4) 

-<o  M 

where  Fj(FK)  is  the  Froude-Krylov  part  of  the  surge  wave-exciting  force. 

Although  simplifing  assumptions  were  made  to  arrive  at  Equation  (3),  their  validity  may 
not  be  justified.  Model  experiments  in  stern-quartering  and  following  waves  have  indicated 
that  large  surge  motion  in  the  order  of  the  wave  orbital  motion  can  be  induced  simultaneously 
with  large  trim.  It  is  quite  plausible  that  coupling  between  the  surge  and  pitch  modes  is  not 
negligible  when  a SWATH  ship  undergoes  large  surge  motion  since  the  vertical  center  of  gravity 
is  located  far  above  the  main  hull  axis.  This  need  to  develop  a more  consistent  hydrodynamic 
theory  for  SWATH  ships  in  following  waves  is  well  recognized,  but  no  such  attempt  will  be 
made  in  this  report. 

The  symmetry  of  the  hull  with  respect  to  the  longitudinal  centerplane  of  a twin  hull 
leads  to  decoupling  of  the  vertical-plane  modes  from  the  horizontal-plane  modes.  Thus,  the 
equations  of  motion  can  be  divided  into  the  following  two  groups: 

Heave  and  pitch  equations; 

(M  + A33)  £3  + B33{3  + C33£3  + A35£5  + B35£5  + C35£5  = F3(e)  e-^' 

, , (5) 

(l5  + A5  5)  £*5  + BS  5^5  + C5S^S  + A5  3^3  + BS  3^3  + C53^3  = F5^  e”^Wt 

where  1^  (=  M^)  is  the  mass  moment  of  inertia  about  the  y-axis. 
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Sway,  roll,  and  yaw  equations; 

+ + B22^2  + ^A24  ~ ^z0)  ^4  + B24^4  + A26^6  + B26^6  = F2^  e 

<>4  + A44>  i4  + B44$4  + C44f4  + <A42  ' Mzo>  *2  + B42*2  + A46*6  + B46*6  = F ^ 

(l6  + A66*  ^6  + B66^6  + A62^2  + B62^2  + A64^4  + B64^4  = F6(c * e_ia,t  (6) 

where  I4  and  I6  are  respectively  the  mass  moment  of  inertia  about  the  x-  and  z-axes,  and  zQ 
is  the  /.-coordinate  of  the  center  of  gravity  of  the  ship. 

The  restoring  coefficients  contributed  by  the  inherent  buoyancy  effect  are  easily 
obtained  by 

C33  = Pg  Aw 
C35  =C53  = PgMw 

_ <7) 
C44  =PgIw4-MgBG 

C55  aP8lw5  -MgBG 


where  p = density  of  water 

g = gravitational  acceleration 
Aw  = waterPlane  area 

Mw  = waterplane  area  moment  about  the  y-axis 
lw4  and  lw3  = moment  of  inertia  of  the  waterplane  area  about  the  x-axis  and  the  y-axis, 
respectively 

BG  = vertical  distance  between  the  center  of  buoyancy  and  the  center  of  gravity 
when  the  ship  is  at  its  mean  position 

As  will  be  seen  later,  there  are  additional  contributions  to  the  restoring  coefficients  from  the 
stabilizing  fins. 

Since  Equations  (5)  and  (6)  are  linear,  the  solutions  of  £k  are  expected  to  be  harmonic 
functions  of  time;  hence  the  equations  can  be  expressed  as  complex  algebraic  equations; 

AlX,=Bi 

A2X2 “ B2 

where 
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ho 

II 

xT 

*3o 

X2  = 

ho 

*50 

^60 

Bl  = 


f3(o 

F5(e> 


B2  = 


F^* 

pU) 

F6(e) 


A.= 


-co2  (M  + A33)  + C33  - j co B 
" A53  + C53  ~Ju’B53 


>33 


-to2  A35  +C35  - jcoB35 
-co2  (I$  + A55)  + C55  - jcoBss 


-co2  (M  + A22)- jcoB22  -CO2  (A24  - Mz0)-jcoB24  -co2  A26  - jcoB26 

-co2(A42-MzQ)-jcoB42  -co2  (l4  + A44)  + C44  - jcoB44  -co2  A46  - jco  B46 

'64  _Jwo64  '*6  ’ A66 

An  inversion  of  the  matrices  provides  the  amplitudes  of  the  motion  by 


-co- A62  - j wB62  -co2  A^4  - jco  BA4  -co'-tl^  + AAft)  - jcoB 


Ik  = l*kol-IStc2+«k,2l''2 

and  the  phase  angles  with  respect  to  the  wave  crest  above  the  coordinate  origin  by 

ak  = tan-'  [-$ks/*kcl 


66 


(8) 


(9) 


for  k = 2.  3,  4,  5,  6. 

The  velocity  and  acceleration  of  the  motion  are  obtained  by  simply  multiplying  the  complex 
amplitudes  of  the  displacements  by  -jco  and  -co2,  respectively.  The  surge  amplitude  and 
phase  c*|  can  be  obtained  similarly  from  Equation  (4). 

The  complex  amplitudes  of  the  absolute  and  relative  vertical  motion  of  a point  (x,  y,  z) 
on  the  ship  are  given  for  the  absolute  vertical  motion  by 


and  for  the  relative  motion  by 


VA)  =S3o  + y*4o-x*5c 

t (R)  = t (A)  _ f 
«v  'v  5o 
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(12) 


where  f0  is  the  complex  amplitude  of  the  incoming  wave  given  by 

► _ A jKo(xcos0-y  sin0) 

f o " A c 

Here  A is  the  wave  amplitude,  K(1  the  wave  number  defined  by  KQ  = 2rr/\  = to02/g,  X the  wave- 
length, gj()  the  wave  frequency,  and  (3  the  wave  heading  angle  with  respect  to  the  x-axis  in 
the  counterclockwise  direction  from  the  x-axis  (=  0 deg  represents  the  wave  heading  from  the 
stern  to  the  bow),  Strickly  speaking.  Equation  (II)  is  an  approximation  since  we  have 
neglected  deformations  of  the  wave  elevation  along  the  body  caused  by  the  diffraction  of 
waves  by  the  body  surface  and  by  the  waves  generated  by  the  body  motion. 

The  complex  amplitude  of  absolute  transverse  motion  of  a point  on  the  ship  is  given  by 

WA)  =*2o  + x*6o-z*4o 

In  the  foregoing,  the  frequency -response  of  the  various  motion  amplitudes  and  phases 
were  given.  The  complex  amplitudes  of  any  motion  quantity  divided  by  the  amplitude  of  the 
waves  are  often  called  transfer  functions  or  frequency-response  functions.  In  principle,  the 
time  history  of  the  system  response  to  any  random  signal  can  be  obtained  for  a linear  system 
by  a convolution  integral  of  the  product  of  the  signal  and  the  inverse  Fourier  transform  of  the 
transfer  function.  A correct  approach  to  ship  motion  in  the  time  domain,  however,  is  not 
so  straightforward  as  described  in  the  foregoing;  this  has  been  pointed  out  by  Cummins.8 
As  is  well  recognized,  a unique  representation  of  the  time  history  of  sea  waves  is  impossible; 
hence  representation  of  sea  waves  has  been  made  through  energy  spectra  from  which  various 
statistical  averages  of  the  wave  conditions  can  be  obtained. 

St.  Denis  and  Pierson9  were  the  first  to  introduce  an  application  of  sea  energy  spectra  in 
conjunction  with  the  transfer  function  to  obtain  various  statistical  averages  of  ship  responses. 
Since  then,  statistical  averages  have  been  used  almost  as  a standard  tool  for  the  investigation 
of  ship  motion  in  irregular  waves.  The  major  underlying  assumptions  for  the  concept  are  that 
the  relationship  between  the  wave  excitation  and  ship  response  is  linear,  that  wave  and  ship 
motion  are  stationary  and  normal  random  processes  with  zero  mean,  and  that  the  spectral 
density  functions  of  waves  and  ship  motion  are  narrow  banded.  If  a sea-energy  spectral 
density  function  which  has  the  dimensional  units  of  |L2T|  is  denoted  by  S(to0),  the  variance 
of  a motion  quantity,  say,  £k.  can  be  obtained  by 

8Cummins,  W.E.,  “The  Impulse  Response  Function  and  Ship  Motions.”  Schiffstcchnik.  Vol.  9,  pp.  101-109 
(1962);  reprinted  as  DTMB  Report  1661. 

9St.  Denis.  M.  and  W.J.  Pierson,  “On  the  Motion  of  Ships  in  Confused  Seas,”  Trans.  SNAME,  Vol.  61. 
pp.  280-357  (1953). 
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As  can  he  seen  in  the  next  section,  the  wave-exciting  forces  k = I,  ...  6,  are  linearly 

proportional  to  the  incoming  wave  amplitude;  hence  if  the  wave  amplitude  is  taken  as  a unit 
value,  the  resulting  motion  amplitudes  obtained  from  the  solutions  of  the  equations  of  motion 
are  already  factored  by  the  wave  amplitude.  In  the  following,  all  transfer  functions  are  under- 
stood to  be  normalized  by  the  wave  amplitude  unless  otherwise  specified.  Thus,  Equation  (13) 
will  be  written  as 


E=/[k)Z  S(wo)dwo 


(14) 


where  (£k  )*  is  often  called  the  response  amplitude  operator  (RAO). 

Of  the  few  mathematical  expressions  for  sea  spectra,  the  most  frequently  used  by  ship 
motion  investigators  are  the  so-called  Pierson-Moskowitz  spectrum10  and  the  Bretschneider 
spectrum.1 1 More  recently  the  trend  is  to  use  actually  measured  sea  spectra  in  order  to 
examine  the  motion  of  a ship  in  a wide  variety  of  sea  conditions.1 2 For  example.  Miles13 
has  provided  a stratified  sample  of  323  sea  spectra  based  on  measurements  obtained  at  Station 
India  in  the  North  Atlantic.  Proper  weighting  factors  for  the  frequency  of  occurrence  for  each 
of  the  spectra  given  in  Miles13  can  be  determined  from  the  useful  wave  statistics  compiled  by 
Hogben  and  Lumb.14  It  follows,  then,  that  the  probability  can  be  determined  that  a ship 
operating  in  the  North  Atlantic  will  exceed  certain  seaworthiness  characteristics,  e.g.,  vertical 
acceleration,  slamming  and  deck  wetness  per  hour,  etc.  Since  these  predictions  are  made 
under  the  assumption  of  the  Rayleigh  probability  distribution  function,  it  is  tacitly  assumed 
that  distribution  of  the  motion  amplitudes  follows  the  Rayleigh  function  and  this  may  not  be 
quite  true  in  some  cases. 


’"Pierson,  W.J.  and  L.  Moskowitz,  “A  Proposed  Spectral  Form  for  Fully  Developed  Wind  Seas.  Based  on  the 
Similarity  Theory  of  S.A.  Kitaigarodskii,”  J.  Geophys.  Res.,  Vol.  69,  No.  24,  pp.  5181-5190  (1964). 

I!  Bretschneider,  C.L.,  “Wave  Variability  and  Wave  Spectra  for  Wind-Generated  Gravity  Waves,”  Beach 
Erosion  Board,  U.S.  Army  Corps  of  Engineers  TM  1 18  ( 1959). 

1 “Hadler,  J.B.  et  al„  “Ocean  Catamaran  Seakeeping  Design,  Based  on  the  Experience  of  USNS  HAYES,” 
Trans.  SNAME.  Vol.  82,  pp.  126-161  (1974). 

1 ’Miles,  M.,  “Wave  Spectra  Estimated  from  a Stratified  Sample  of  323  North  Atlantic  Wave  Records,” 
National  Research  Council,  Division  of  Mechanical  Engineering  Report  LTR-SH-118  (1971). 

l4Hogben,  N.  and  F.E.  Lumb,  “Ocean  Wave  Statistics,”  H.M.  Stationery  Office,  London  (1967). 


10 


The  variance  for  the  velocity  and  the  acceleration  of  the  kth  mode  motion  can  be  easily 


obtained  bv 


Ev=  / (t*>4k)2  S(w0)dw0 


Ea  =y  <w2  4k)2  S(w0)dw0 


Both  the  Pierson-Moskowitz  spectrum  and  the  Bretschneider  spectrum  can  be  expressed 


in  the  form 


S(w0)  = — ^ exp  (-C2/ou04) 


0.78  - Pierson-Moskowitz 


I 48  7. 06  H2 


- Bretschneider 


in  which  Hs  is  the  significant  waveheight  in  metres,  T0  is  the  modal  period  in  seconds,  and 

3. 12/HS2  -Pierson-Moskowitz 

C2  = 

1948. 24/T  2 - Bretschneider 


In  the  application  of  these  empirical  formulas,  caution  is  necessary  in  taking  the  integral 
limit  to  infinity  since,  as  can  be  seen  from  Equation  (16),  cj4  S(ojo)~*  1/ojo  as  hence 

the  integral  behaves  like  a logarithmic  function  with  an  infinite  argument  e.g.,  since  the  RAO 
of  relative  vertical  motion  approaches  unity  as  wQ  goes  to  infinity,  the  spectrum  of  the 
acceleration  of  relative  motion  at  high  frequencies  would  behave  like  l/co0.  Since 
Equation  (16)  is  an  empirical  formula  based  on  the  recorded  ocean  waves  which  usually  yield 
uncertain  data  at  high  frequencies  because  of  the  filtering  process,  the  high-frequency  end  of 
the  formula  is  not  reliable  as  Pierson  points  out  (see  pages  86-89  of  Reference  15).  Thus,  for 
the  prediction  of  motion  of  ships  with  lengths  greater  than  30  m.  the  high  frequency  limit  in 
the  integration  should  be  replaced  by  something  like  3.0  rad/'sec  which  corresponds  to  about 
a 6-m  gravity  wavelength  in  deep  water. 

15Piersor,,  W.J.,  “The  Theory  and  Applications  of  Ocean  Wave  Measuring  Systems  at  and  below  the  Sea  Surfaces, 
on  the  Land  from  Aircraft,  and  from  Space  Craft,”  NASA  Contractors  Report  NASA  CR-2646  (1976). 
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Various  often  used  statistical  averages  can  be  expressed  in  the  form 


Average  Amplitude  = (17) 

Here  E0  is  the  variance  of  a particular  motion  amplitude.  C=  1.253  provides  the  average, 

C = 2.0  provides  the  one-third  highest  average  or  significant  average  and  C = 2.546  provides 
the  one-tenth  highest  average. 

From  the  Rayleigh  law  of  probability  distribution,  the  probable  number  of  slams 
sustained  per  n sec  by  the  main  hull  bottom  or  the  cross-deck  bottom  of  a SWATH  ship  at  a 
given  location  can  be  given  by 


Ns  = T* 
s 2j t 


2E(R)  2Ev(R) 


The  superscript  (R)  denotes  the  relative  motion  of  the  location  of  interest,  CQ  is  the  vertical 
distance  from  the  calm  waterline  to  either  the  main  hull  bottom  or  the  cross-deck  bottom  in 
the  same  dimensional  unit  used  for>/ E(R)  , and  VT  is  the  threshold  velocity  that  incites 
slamming.  The  value  of  VT  can  differ  from  case  to  case  and  should  be  given  in  the  same 
dimensional  unit  as  \/ . In  many  cases  the  values  of  VT  are  unknown,  and  if  VT  is  set 
to  zero.  Equation  (18)  then  gives  either  the  number  of  hull  bottom  emergences  or  water 
contacts  of  the  cross-deck  bottom  per  n sec.  If  we  set  VT  = 0 and  take  CQ  to  be  the  deck 
height,  then  Equation  (18)  provides  the  probable  number  of  occurrences  of  deck  wetness  per 


For  a design  or  operational  criterion,  it  is  also  of  interest  to  know  the  probable  extreme 
value  a ship  may  encounter  in  given  sea  environments.  Ochi16  has  shown  that  the  extreme 
value  in  amplitude  expected  in  n observations  yn  can  be  expressed  by 


for  small  6 

where  T is  the  time  in  hours  during  which  an  extreme  sea  environment  may  persist,  and  6 is 
given  a value  of  0.01  if  the  design  goal  calls  for  a 99-percent  assurance  that  the  extreme 
amplitude  yn  will  not  be  exceeded.  For  5 = 1. 7n  represents  the  “most  probable  extreme 
value”  in  amplitude.  For  a large  number  of  observations,  the  probability  that  the  extreme 
value  will  exceed  7n  is  63.2  percent. 

l6Ochi,  M.K.,  “On  Prediction  of  Extreme  Values,”  J.  Ship  Res.,  Vol.  17,  No.  1,  pp.  29-37  (1973). 
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EVALUATION  OF  HYDRODYNAMIC  COEFFICIENTS 

The  hydrodynamic  coefficients  in  the  equations  of  motion  will  be  divided  into  three 
groups.  The  first  consists  of  those  coefficients  which  can  be  obtained  under  the  assumption 
of  potential  flow,  the  second  consists  of  those  coefficients  mainly  associated  with  the  viscous 
nature  of  the  fluid,  and  the  third  is  associated  with  the  hydrodynamic  lift  generated  by  the 
stationary  fins. 

The  derivation  of  these  coefficients  involves  application  of  various  assumptions  and 
approximations,  and  rigorous  justifications  are  lacking  for  some  of  them.  This  section  attempts 
to  describe  as  well  as  possible  the  underlying  assumptions  involved  in  deriving  each  coefficient. 
There  are  still  unsolved  problems  to  be  looked  into  before  more  satisfactory  justifications  can 
be  offered  for  deriving  some  of  the  coefficients,  especially  those  in  the  second  and  third  groups. 

COEFFICIENTS  OBTAINED  UNDER  POTENTIAL  FLOW  ASSUMPTION 

The  added  mass  coefficients  Ay,  the  damping  coefficients  contributed  by  the  motion- 
generated outgoing  waves  By,  and  the  wave-exciting  forces  F*e)  for  i.  j = 1 . 2.  3,  4.  5,  6 belong 
to  this  group.  The  basic  solution  is  obtained  under  an  assumption  of  two-dimensional 
potential  flow  at  each  cross  section  of  the  ship.  The  boundary-value  problem  is  solved  for  the 
velocity  potential  functions  for  infinitely  long,  semisubmerged,  horizontal  twin  cylinders 
having  cross  sections  identical  to  the  cross  section  of  a SWATH  undergoing  heave,  sway,  or 
roll  oscillation.  This  is  done  by  the  method  of  source  distribution  which  is  described  in  detail 
by  Lee  et  a I.1 7 for  heave  oscillation;  an  extension  of  the  method  was  made  later  for  sway  and 
roll  oscillations.  The  validity  of  the  two-dimensional  solution  was  checked  and  comparisons 
between  the  theoretical  and  experimental  results  presented17  for  rectangular,  circular,  and 
triangular  twin  cylinders. 

Application  of  the  source  distribution  method  to  oscillating  twin  cylinders  results  in  two 
distinct  singular  solutions  at  certain  discrete  frequencies  of  oscillation.  One  stems  from  the 
mathematical  failure  associated  with  the  solution  of  the  Fredholm-type  integral  equations  in 
the  course  of  determining  the  source  strengths,18  and  the  other  stems  from  the  physical  reality 
that  standing  waves  are  trapped  between  the  two  cylinders  at  certain  frequencies  of  oscillation 
given  approximately  by 

l7Lee,  C.M.  et  al.,  “Added  Mass  and  Damping  Coefficients  of  Heaving  Twin  Cylinders  in  a Free  Surface,” 
NSRDC  Report  3695  (1971). 

i8John,  F.,  “On  the  Motion  of  Floating  Bodies:  II.  Simple  Harmonic  Motions,”  Commun.  Pure  Appl.  Math., 
Vol.  13,  pp.  45-101  (1950). 
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u>  = 


' gnf 
(b/a  -1 1 


for  n = 1,2,.. 


where  b is  one-half  the  distance  between  the  centerline  of  each  cylinder  and  a is  the  half-beam 
of  the  individual  cylinder  at  the  waterline. 

The  former  type  of  singular  solutions  can  be  removed  by  extending  the  source  distribu- 
tion onto  the  waterline  inside  the  cylinders  and  imposing  either  wall  condition  or  the  condi- 
tion of  vanishing  velocity  potential  on  the  inner  waterline.  The  validity  of  this  method  can 
be  demonstrated  b>  numerical  results6, 19  and  yet  no  vigorous  mathematical  proof  has  been 
established  for  the  existence  and  the  uniqueness  of  the  solution  obtainable  by  this  method. 

I here  appears  to  be  no  way  to  remove  the  latter  type  of  singular  solutions  except  by  a full 
three-dimensional  solution.  In  practice,  the  range  of  frequencies  of  interest  for  motions  of 
ships  of  lengths  gnuter  than  200  ft  lies  below  the  frequencies  at  which  the  singular  behavior 
of  the  two-dimensional  solution  occurs.  However,  caution  is  necessary  when  the  loading  on  a 
ship  due  to  waves  is  computed,  since  the  peak  loading  may  occur  in  a higher  frequency  range. 

To  obtain  the  three-dimensional  coefficients,  the  sectional  hydrodynamic  coefficients  of 
the  two-dimensional  solution  are  integrated  lengthwise.  For  the  case  of  a ship  without 
forward  speed,  the  strip  approximation  described  in  the  foregoing  may  still  be  an  acceptable 
approximation  in  a practical  sense;  however,  it  appears  that  for  a ship  with  forward  speed, 
the  lengthwise  flow  disturbances  generated  by  the  forward  speed  may  immediately  invalidate 
the  two-dimensional  flow  assumption  at  each  cross  section.  The  long  controversy  on  this 
particular  point  still  continues. 

Ogilvie  and  Tuck20  presented  a more  consistent  and  rational  theory  for  forward  speed 
effects  on  the  hydrodynamic  coefficients  which  are  derived  on  the  basis  of  strip  theory. 

Faltinsen21  demonstrated  a better  correlation  of  the  theoretical  results  obtained  from 
the  Ogilive  and  Tuck  theory20  with  experimental  results  for  a few  of  the  hydrodynamic 
coefficients  which  can  be  strongly  influenced  by  forward  speed  of  a ship.  There  is  no  doubt 
that  computations  of  these  coefficients  according  to  Ogilvie  and  Tuck  are  much  more  tedious 
and  expensive  than  computations  according  to  conventional  strip  theory.5  Yet,  conclusive 


1 ^Ohmatsu.  S.,  “On  the  Irregular  Frequencies  in  the  Theory  of  Oscillating  Bodies  in  a Free  Surface,”  Ship 
Research  Institute  of  Japan  Report  48  (1975). 

2nOgilvie,  T.F.  and  F..O.  Tuck,  “Rational  Strip  Theory  of  Ship  Motions;  Part  1,"  University  of  Michigan, 
College  of  Engineering  Report  Old  (1969). 

21  Faltinsen,  O.,  “Numerical  Investigation  of  the  Ogilivie-Tuck  Formulas  for  Added  Mass  and  Damping 
Coefficients,”  J.  Ship  Res.,  Vol.  18,  pp.  73-84  (1974). 
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evidence  does  not  seem  to  exist  that  this  new  theory  will  necessarily  improve  the  prediction 
of  ship  motion.* 

Although,  as  pointed  out  by  Ogilvie  and  Tuck,  some  argument  may  be  made  regarding 
the  consistency  of  the  perturbation  expansion,  conventional  strip  theory  will  be  used  in  this 
report.  Details  of  the  development  of  the  theory  are  given  in  Appendix  A and  only  the  final 
expressions  obtained  there  will  be  given  here. 

Let  the  velocity  potential  function  1>(x,  y,  z,  t)  whose  gradient  represents  the  velocity 
field  in  the  fluid  region  disturbed  by  incoming  waves  and  motion  of  a ship  be  expressed  in 
the  form 

6 

<l>  = -U  x + 0s(x,  y.  z)  + (0|(x,  y,  z)  + 0D  + £jo  0j)  e~iwt 

i=  I 

Here  0s  is  the  steady  wave-resistance  potential,  0(  and  0q  are  respectively  the  complex 
potentials  representing  the  incoming  and  diffracted  waves.  0j  is  the  complex  potential  repre- 
senting the  fluid  disturbance  by  the  ship  motion  in  the  it h mode,  and  £jo  is  the  complex 
amplitude  of  the  displacement  in  the  it h mode.  Then,  the  hydrodynamic  coefficients  can  be 
given  by 


Aik  = Rej  ~ /dx  [ *,:N 

L JL  ~/C(x ) 


2U  , 2U 

;N(y.  z:x)  + — 03 N 6i5  - — 0jN  5i6  } <t> 


kd« 


(20) 


Bik  = lr^ 


~u  ^dxJ^xJ0iN(y’Z:X)  + j^03N6i5-j^^2N6i6j^kd8  <2,) 


for  i.  j = 2.  3,  4,  5,  6 (sec  Equations  (76)  and  (77)  in  Appendix  A). 


In  Equations  (20)  and  (21),  Re^  and  Iirij  mean  the  real  and  the  imaginary  part  of  what 


follows 


: Jl  iS 


the  lengthwise  integral; 


■J 


is  the  integral  along  the  immersed  contour  of 


C(x ) 


the  section  located  at  x when  the  ship  is  at  its  mean  position;  6jk  is  the  Kronecker  delta 
function;  <t>[  is  the  two-dimensional  representation  of  0;  (strip  assumption);  and  subscript  N 
denotes  the  normal  derivative  in  the  y-z  plane,  the  positive  sense  of  the  two-dimensional  unit 
normal  vector  N being  pointing  into  the  body.  For  later  use.  the  sectional  added  mass  a^lx) 
and  the  wavemaking  damping  b^tx)  are  defined  by 


•Since  as  their  first  effort.  Ogilvie  and  Tuck  covered  only  thp  coefficients  on  the  left-hand  side  of  the 
equations  of  motions  (see  Equation  (5)),  motion  computation  cannot  be  performed  unless  a consistent 
theory  is  also  developed  to  cover  the  wave-exciting  terms. 
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bjj(x)  = lm- 


(22) 


(23) 


for  i = 2,  3 and  4. 

The  wave-exciting  terms  are  expressed  by 

F,(e)  = jcJ0jO  f dx  f n,  0,  df 
JL  -'C(x) 


(24) 


Ffe)  = p y"dx  J £u0  Nj+  ( 03  6i5-~^2  6ib)  alTl  ^1  dC  (25) 

L C ( x ) 

for  i = 2,  3,  4,  5.  6.  Here  nj  is  the  x-component  of  the  unit  normal  vector  on  the  ship  surface 
pointing  into  the  body,  Nj  is  the  ith  component  of  the  two-dimensional  unit  normal  vector 
such  that  N-,  and  N3  are  the  y-  and  z-components, 

N4  = y N,  - z N2 

N5  =-xN3 

N6  = x N2 


and  0 j is  given  by 


and 


0I  = -:^exp(Koz+jKox  cos  (3  - j K0y  sin  0) 


0,N  = K0  (-jN2  sin  0 + N3)  0| 


(26) 


in  which  the  notations  are  as  defined  under  Equation  (12). 

If  the  hull  geometry  is  given  by  z = h(y;x)  the  definition  of  the  components  of  the  unit 
normal  vector  can  be  given  as 

1 ( dy(x)  dz(x)\  * 


ni  ‘ FI  ( 


y dx  dx  / 


•For  the  main  hull  of  a SWATH  ship  having  the  form  of  body  of  revolution,  we  can  obtain  n,  by 
R'(x)/(1  + R'2)1/2  where  R(x)  is  the  radius  of  the  cross  section  of  one  hull  and  R'(x)  = dR/dx. 
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where  N = y'  1 + h “ . If  the  offsets  of  the  hull  are  given  at  several  waterlines  at  each  station 
of  a ship,  the  hull  derivatives  dy/dx  and  dz/dx  may  be  obtained  numerically  by  use  of  the 
three-point  Lagrangian  interpolation  rule  along  the  same  waterline.  As  can  be  seen  in  Figure 
2.  if  we  know  the  tangent  angle  a at  a given  point  of  a cross  section  contour,  then  at  a given 
point  of  a cross  section  contour,  we  have 


hy  = tan  a 


Hence 


. N’t  = ~sin  a and  N,  = cos  a 

cos  q -> 
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Table  I shows  the  expressions  for  the  hydrodynamic  coefficients  given  in  terms  of 
Equations  (20)  through  (2b).  Note  that  the  forward  speed  effect  on  the  hydrodynamic 
coefficients  appears  distinctly  as  multiplication  factors  to  some  coefficients.  However,  the 
sectional  hydrodynamic  quantities  a^  and  for  i = 2,  3,  and  4 are  functions  of  the  geometry 
of  the  cross  section  as  well  as  the  encounter  frequency,  and  since  the  encounter  frequency  is 
a function  of  the  forward  speed,  the  forward-speed  effects  are  also  implicitly  imbedded  in 
other  hydrodynamic  coefficients. 

To  facilitate  reader’s  undeistanding,  the  derivation  of  A$3  and  B66  will  be  illustrated  in 
the  following.  From  Equation  (20) 


From  Equation  (68)  in  Appendix  A, 


N 


3 


03N 

jw 


05fyj  j oj  Nj  + UN3 


From  the  slender-body  assumption. 


Hence 


and 


N_s  = - x N3 

«5N  =j“(x  +j^)  N3=~  (X+jS)  03N 
05=-(x+jz;)*3 


Substitution  of  this  relation  into  the  expression  of  A53  above  yields 
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TABLE  I - HYDRODYNAMIC  COEFFICIENTS  OBTAINED 
UNDER  POTENTIAL-FLOW  ASSUMPTION 


A22  = / a22 (x) dx 

JL 

B b2^  dx 


^24  - J" a24  dx 


®24  b24  dx 


A35  = -/xa33  dx-^B33 

CO^ 


-I 


35  =-/  x b33  dx  +UA33 


A44  = / a44  dx 


B44  “ / b44  dx 


A26  J x a22  dx  + ~ 5 B22 

u> 


^26  7"" 
>26=y~> 


J x b-^  2 dx  U ^22 


A42  A24  ^7 a24  dx 

B42  = B24  = /b24  dx 


A33  / a33  dx 


A46  = / x a24  dx  +“^  B24 

CO" 


B33  ~ I b33  dx 


B46  " / x b24  dx  " U A42 


r , uz 

A55=  / x a33dx+“  A33 
-/  co 


B55  =/x2  b33dx+— 2 B33 

^ CO 


A53  " x a33  dx  + ~ B33 
'/  ur 


B53  ~ J x b33  dx  " U A33 


-P 


r ■>  u2 

A66  ~ j x a22dx+~  A22 

^ CO 


f , u2 

~ J x“  dx  H - B-)2 

co"  ~ 

A62  = / x a22  dx  2 B22 
co" 


B62  ~ J x B'>2  dx  + U A' 


A64  J x a24  dx  2 B24 
co" 

B64  =^J~ x b24  dx  + ^ A24 
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TABLE  1 (Continued) 


F,(e>  = 

PgA  / 
JL 

i k x cos  p r 

dx  e nl  / (~j  ^2  sin  i 

-'C(x) 

8 + NJ,el<"<!-iyim'!,d» 

F\(e>  = 

-jpg  A 

CO 

0 * 

/dxejK»’ 

'L 

t cos  0 r | 

Ax)l 

jw0  N2  + K0(-j  N2  sin  (3  + N3  ) 02 

F3<e>  = 

-j  PgA 

CO 

0 * 

/dxejKo 

a 

x cos  (3  r 

Ax) 

{ J^o  N3  + 

KQ  (-j  N2  sin  P + N3)  03 

eK„(z-jy™|3ld|! 

F4,e>  = 

-jpg  A 

CO 

0 ' 

/dXeiK» 

x cos  (3  r 

Ax) 

jjco0(yN3 

- z N2 ) + KQ  (-j  N2  sin  /: 

! 

+ N3 ) 04  J 

eKo(Z-jysin% 

F5<e>  = 

-jpgA 

CO 

o % 

C iK  ) 

/ dx  eJ  ° 

a 

i cos  (3  r | 

Ax)  1 

J-jWoxN3 

+ KQ  (-j  N2  sin  P + N3) 

.1  K (z  -jy  sin  (3) 

\jw  ) 

p3  | e ux 

f6w. 

-jpgA 

CO 

0 ~ 

r j k > 

/ dx  e ° 

'l 

c cos  (3  A"  1 

Ax)' 

J jco0  X N2  ■ 

- KQ  (-j  N2  sin  P + N3 ) 

* (jw"X)  ^ eK°<Z"iySin^ 


0- >>  e 


dC 
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The  definitions  given  by  Equations  (22)  and  (23)  can  be  used  to  write 


(x ) dx 


Procedures  similar  to  the  foregoing  together  with  the  relations 


N2  = ~02N^<*O>  N6  = x N2,  and  <t>'6  = 


can  be  used  to  derive 


VISCOUS  DAMPING  COEFFICIENTS 

Hydrodynamic  coefficients  derived  under  the  potential-flow  assumption  (as  shown  in  the 
preceding  section)  have  been  found  to  provide  satisfactory  motion  predictions  for  most 
conventional-type  ships  wherein  wavemaking  damping  predominated.  However,  as  in  the  case 
of  predicting  roll  motion  of  a surface  ship,  when  wavemaking  damping  is  no  longer  the 
dominant  contributing  factor  to  the  overall  damping,  the  damping  contributed  by  the  viscous 
effects  of  the  fluid  has  to  be  taken  into  account. 
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The  scmisubmersible  SWATH  configuration  does  not  generate  large  surface  waves  when 
it  oscillates  in  the  vertical-plane  modes.  This  means  that  the  wavemaking  damping  of  SWATH 
ships  in  the  vertical-plane  modes  (i.e.,  heave,  pitch,  and  surge)  is  relatively  small  compared 
with  that  of  conventional  ships.  When  the  viscous  effects  contributing  to  the  damping  are 
neglected,  the  computed  motion  amplitude  in  the  neighborhood  of  resonance  is  similar  to 
that  of  a typical  underdamped  linear  system,  i.e.,  a narrowly  tuned,  high  spiked  motion  at 
the  resonant  frequency.  Figure  3 shows  such  behavior  for  a SWATH  configuration  proceeding 
in  regular  head  waves  when  computed  with  the  potential-flow  hydrodynamic  coefficients 
alone.  Model  Experimental  results  are  included  in  the  figure  to  indicate  the  errors  caused 
when  the  predictions  underestimate  damping  coefficients.  The  need  for  a more  reasonable 
evaluation  of  damping  coefficients  is  obvious;  however,  discouragingly,  there  is  very  little  one 
can  do  to  improve  the  evaluation  theoretically. 

Encouraged  by  the  remarks  made  in  Pien  and  Lee,6  Kim22  tried  to  include  the  forward- 
speed  effect  more  rigorously,  and  used  a thin-ship  approach  in  his  evaluation  of  damping 
coefficients.  Figure  4 compares  his  results  with  those  obtained  by  the  strip  theory  (described 
in  the  preceding  section)  and  by  experimental  results  for  a SWATH  model  designated  as 
SWATH  I (previously  called  Modcat  1).  Kim’s  results  do  not  differ  significantly  from  the 
experimental  results,  yet  when  the  motion  based  on  these  damping  coefficients  was  computed, 
the  unrealistically  large  motion  that  characterized  Figure  3 were  still  present.  This  implies 
that  even  experimentally  obtained  damping  coefficients  are  not  sufficient  to  yield  reasonable 
motion  predictions.  Thus  we  are  led  to  conclude  that  forced-oscillation  experiments  of  a 
model  with  small  amplitude  of  motion  do  not  yield  realistic  values  of  damping  for  a model 
that  undergoes  large  motion  near  reasonance. 

If  the  foregoing  conclusion  is  correct,  the  damping  must  be  nonlinear  with  respect  to  the 
motion  amplitude  when  it  exceeds  a certain  magnitude.*  Figure  5 presents  unpublished 
experimental  results  obtained  at  DTNSRDC  by  Wahab  by  vertically  oscillating  a two-dimen- 
sional cylinder  with  a bulbous  bow  cross  section  at  zero  speed.  The  theoretical  results  shown 
in  the  figure  were  obtained  by  using  the  source  distribution  method  developed  by  Frank  for 
single  cylindrical  forms.23  Here,  a large  discrepancy  between  the  potential-flow  theory  and 
experiment  is  obvious,  unlike  the  case  for  a three-dimensional  body  shown  in  Figure  4. 

22Kim.  Ki-Han.  “Determination  of  Damping  Coefficients  of  SWATH  Catamaran  Using  Thin  Ship  Theory.” 
Mass.  Inst.  Technol.  Dept.  Ocean  Eng.  Report  75-4  (1975). 

*A  linearity  check  for  the  damping  coefficient  with  respect  to  the  motion  amplitudes  was  made  in  the 
experiment  as  shown  in  Figure  4. 

23Frank,  W..  “Oscillation  of  Cylinders  in  or  below  the  Free  Surface  of  Deep  Fluids,”  NSRDC  Report  2375 
(1967). 


HEAVE 


Figure  3 — Heave  and  Pitch  Motion  of  SWATH  4 in  Regular  Head  Waves  at 
20  Knots  Without  Viscous  Damping 


A forced  heave -oscillating  experiment  with  a SWATH  model  (some  results  of  which  are 
shown  in  Figure  4)  was  conducted  at  DTNSRDC  by  Stahl  (as  yet  unpublished)  under  various 

conditions.  The  forced-vertical  oscillation  test  were  made  with  the  models  of  a single  hull  of 
two  different  scale  ratios  and  twin  hulls,  with  three  to  four  different  amplitudes  of  oscillations 
to  within  0.16  of  the  draft,  four  Froude  numbers  (Fn  = 0,  0.2,  0.4  and  0.6),  and  a wide  range 
of  frequencies.  The  nonlinear  effect  on  the  damping  coefficient  was  not  significant  except 
perhaps  for  low  frequencies  at  high  speeds.  The  scale  effect  for  a single  hull  appeared  signifi- 
cant only  at  the  high  frequencies  which  are  outside  the  range  of  interest  for  SWATH  ships. 
Speed  effect  does  not  appear  significant  except  in  the  very  low  frequency  range.  The  mutual 
blockage  effect  of  twin  hulls  on  damping  appeared  to  be  insignificant.  This  experiment,  in 
fact,  has  demonstrated  that  the  mathematical  model  derived  under  the  assumption  of  small 
motion  amplitude  can  be  justified  for  motion  amplitudes  within  the  range  in  which  the 
experiment  was  conducted. 

It  has  frequently  been  observed  that  SWATH  models  have  rapid  transient  decays  when 
they  have  forward  speeds.  Thus,  one  might  be  led  to  conclude  that  a SWATH  configuration 
possesses  a critical  damping  when  at  speed.  However,  SWATH  models  at  speed  have  exhibited 
a resonant  motion  when  they  are  subjected  to  continuous  harmonic  wave  excitations. 

The  foregoing  description  is  intended  merely  to  emphasize  the  complexity  of  the  problem 
we  are  attempting  to  solve  and  thus  demonstrates  that  a simplified  theory  could  hardly  be 
expected  to  provide  more  than  qualitative  agreement.  Left  without  an  alternative,  it  was 
decided  to  follow  an  empirical  approach  to  determine  the  supplemental  damping  required  to 
predict  the  motion. 

According  to  Thwaites,24  experimental  results  of  side  forces  generated  on  slender  bodies 
with  moderate  trim  angle  showed  that  the  vertical  force  Fy  can  be  expressed  in  the  form 

1 s 

Fv  = — p U“  Ap  sin  a | sin  a | (a0  I cot  a|  + CD)  (27) 

Here  Ap  is  the  projected  plane  area  of  the  body  on  a horizontal  plane,  a is  the  trim  angle,  and 
a0  and  CD  are  real  constants:  a0  is  often  called  the  viscous-lift  coefficient  and  CD  the  cross- 
flow  drag  coefficient.  The  value  of  a0  found  from  the  experiments  on  airship  models  with 
circular  or  polygonal  cross  sections  was  about  0.07  while  Cq  seemed  to  lie  beween  0.4  and 
0.7.  An  experiment  by  Allen  and  Perkins25  showed  that  CD  increases  as  the  slenderness 
ratio  increases. 

^Thwaites.  B.  (Editor).  “Incompressible  Aerodynamics,”  Oxford  University  Press,  pp.  405-421  (1%0). 

25Allen,  H.J.  and  E.E.  Perkins,  “A  Study  of  Effects  of  Viscousity  on  Flow  over  Slender  Inclined  Bodies  of 
Revolution,”  NACA  Report  1048  (1951). 


For  a harmonically  oscillating  body  in  the  vertical-plane  modes  in  regular  waves  with  a 
constant  forward  speed,  the  foregoing  expression  will  be  assumed  to  take  the  form 

Fv  = ~ p Ap  (U2  a0  a + CD  w | w | ) 

where  w is  the  relative  fluid  velocity  with  respect  to  the  body  and  given  by 
w=*3-x{5  +yt4-fv(x,  ±b(x),-d,(x»=  U(a-£5) 

I he  angle  of  incidence  ot  flow  at  a cross  section  of  the  body  at  x can  be  expressed  as 

<*  = £5  + (£3  - x £5  + y £4  - fv(x,  ± b(x),  - d,  (x ))/U  (28) 

where  fv  = vertical  velocity  of  the  fluid  induced  by  the  incoming  wave 

b(x)=  transverse  distance  from  the  x-axis  to  the  midpoint  of  the  beam  of  one  hull 
d j ( x ) = depth  to  the  maximum-breadth  point  at  a cross  section 

The  expression  given  above  tacitly  assumes  that  a is  small  and  the  diffraction  of  the  incident 
wave  can  be  neglected.  It  should  be  noted  that  the  term  y£4  will  not  contribute  to  the 
vertical  force  because  of  the  symmetry  of  ship  geometry  with  respect  to  the  centerplane. 

If  the  body  is  undergoing  motion  in  the  horizontal-plane  modes,  the  angle  of  flow 
incidence  at  the  mean  depth  of  the  cross  section  of  x can  be  expressed  as 

T?  = -£6  +(£2  +x*6  +d2(x)^4~fH(x’±b(x)’-d2<x))/U  (29) 

where  d-,(x)  is  one-half  the  draft  at  the  cross  section  at  x and  (x,  y,  z)  is  the  transverse 
velocity  of  the  fluid  induced  by  the  incoming  wave. 

Now,  similar  to  the  potential-flow  case,  the  strip  concept  will  be  introduced  to  Equation 
(27)  together  with  the  assumption  of  small  a and  no  viscous  interactions  between  the  two 
hulls.  The  vertical  force  induced  on  the  twin  hulls  can  then  be  obtained  by 

Fv  =■?  U2  /~Bm(x)  V <a0  a/x)  + Cploij  - £5) - £5  I ) dx 
i = I 

where  Bm(x)  is  the  maximum  beam  of  the  submerged  cross  section  of  one  hull  at  x and  a, 
and  a2  arc  respectively  the  angle  of  incidence  on  the  port  and  starboard  hulls.  With  the 
following  new  notation  defined  by 

*1  =ils+ilp- 

^ls  = $3  ~ x £5  ■ fv  (x,-b(x).-d,(x)) 
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and 


Zip  = $3  - -?v(x*  Njc),  -dj(x)) 


Fv  can  be  written  as 


Fv  =4  /Bm(x)|2a0U2£5+a0UJ.i  + CD  (z)s  | z,s  | + z,  I z,  | ] dx  (30) 
" '/L 

Similarly,  the  horizontal  force  FH  can  be  expressed  as 

Fh  = ^2  /*d<x)^  <at>T7j(x)  + CD(r7i  + €6 » I i?j  + £6  l)dx 

' */L  j=l 

= : ^dI':aou2^6 + aouyi + cd (yis 1 vis 1 + yi P 'yip 


dx  (31) 


where  d(x)  is  the  draft  of  the  cross  section  at  x. 


vi  = yis  + yi. 


yls  = k2  + x^6  +d2(x)^4~^H  <x. -b(x), -d2(x)) 
y|p  = ^2+x^6+d2^4-^H  (X'  b(x)'  -d2(x)) 


in  Equations  (30)  and  (31).  a0  and  CD  are  assumed  to  be  constant  over  the  length  of  the 
ship. 

The  moments  contributed  by  the  flow  incident  angles  can  similarly  be  given  for  roll 
moment  by 


dx 


Mr  =?  / Bm(x)b(x)[a0Uzr  + CD(zrs|zrsl  + zrp  |z  ' 

z Jl. 

+ ^ /'d(x)d^(x)|-2a0U2^6 +a0Uy,  +CD  (yls  I yls  I + y,  I y,  1)1  dx  (32) 

" -'l 


where 


zr  = zrs  + zrp 


zrs  = b(x)j4  +fv(x.-b(x).-d,(x)) 
zrp  = b(x)£4  - fv  (x,  b(x).  — d ! ( x )) 
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<■ 


, 


‘ " “ 1,11  ' ' 


For  pitch  moment,  the  equation  is 

Mp=-|j^xBm(x)12aou2C5+aoUi  + CD(iisU|sl+iip|i1pl)l  dx  (33) 
and  for  yaw  moment: 

MY=|  ^"xd(x)[-2a0U2{6  +aQUy,  + CD  (y, s I y, s I + y , p I y,  p I ) ] dx  (34) 

The  fluid  velocity  induced  by  the  incoming  wave  is  obtained  from  Equation  (26)  by 
90]  . 

?v(x,  y,  z,  t)  = — e'-l00’ 
v 3z 

= e'Jwt  ^ ^-^exp(K0z  + j KQx  cos  0 - j KQ  y sin  0)j 
= -jcu0  A exp(KQz  + j Kqx  cos  0-  j KQ  y sin  0)  e-^*  (35) 

and 

90]  . , 

f„(x,  y,  z,t)=  — e-Jwt 

= -u>0  A exp  (K0z  + j KQx  cos  0- j KQy  sin  0)  e~ju,t  (36) 

The  cross-flow  drag  terms  in  the  foregoing  expressions  are  nonlinear;  hence  they  cannot 
be  directly  introducted  into  the  linear  equations  of  motion.  By  the  rule  of  Fourier,  it  can  be 
shown  that 


cos  9 | cos  9 | = • 


L 

i 


2n 

cos  9 1 cos  9 | cos  n 9 d0 


2tt 

cos2  n 9 69 


cos  n 9 


= Aq  + A]  cos  9 + A-,  cos  26  + . . . 
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where 


An  = 0 for  n even 

A = < — 1 » n for  n odd 

n(n*  - 4)tt 

i -JL  x - JL 

A'  ‘ 3tt • A3  “ 1 5tt 

Thus,  for  any  harmonic  motion  given  by  x = xt)  cos  cat,  an  approximation  of 

8 

x | x | — — co  x x (37) 

STT  ° 


can  be  established.  This  approximation  is  often  called  the  equilinearization  method,  and  it  is 
used  frequently  for  a dynamic  system  with  weakly  nonlinear  behavior. 

Substitution  of  the  cross-flow  drag  terms  in  the  foregoing  equations  for  the  forces  and 
moments  by  the  equilinearized  form  requires  a prior  knowledge  of  the  amplitude  of  motion 
such  as  x0  in  Equation  (37).  In  the  present  work,  the  amplitude  of  motion  is  obtained  by 
solving  the  equations  of  motion  without  the  cross-flow  damping  terms:  the  equations  of 
motion  are  then  resolved  with  the  equilinearized  damping  terms,  and  the  process  is  repeated 
until  a reasonable  convergence  on  the  motion  amplitude  is  obtained. 

The  viscous  damping  coefficients  and  the  wave-exciting  coefficients  can  be  derived  from 
the  foregoing  expressions  of  the  forces  and  moments  generated  by  the  viscous  effects.  These 
coefficients  are  expressed  with  an  asterisk  attached  to  the  notation  and  are  shown  in  Table  2. 
The  definitions  of  the  new  notation  used  in  Table  2 are  given  by 


The  restoring  coefficients  such  as  C^.  and  cannot  be  incorporated  in  Equation  (16) 
since  such  restoring  terms  are  not  present.  Thus,  they  will  be  lumped  with  the  added  inertia 
term  by  equating  C’*6  = - A^/co2.  = - A^/ta2,  and  = - A^/oa2.  As  long  as  the 

equations  of  motion  in  linear  form  are  solved  in  the  frequency  domain,  such  an  exchange  of 
the  terms  between  the  inertial  and  restoring  coefficients  with  the  factor  (-ca2)  will  not  alter 
the  solutions.  However,  this  exchange  of  the  terms  should  not  be  made  in  the  time-domain 
analysis  of  the  motion  and  in  the  stability  analysis  of  the  ship. 
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TABLE  2 - VISCOUS  DAMPING  COEFFICIENTS 


B *2  =Pa°Ujfd(X)dX  +P  ^CD^lo 

b24 = p ^y d ^ dx + p Tv  cd  y o 
B2*6  =Pa0Uyxddx+p^CD  yyl0> 
C2*6  =- A2*6/w2  =-pa0  J' d dx 

B33=PaoU  dx^CD/^ 


(x) d(x ) dx 


d d2  dx 


Bm  dx 


'3*5  = ~PaoU  /KBmdx-f3j;C0  J X'loBmdx 


C,%  =Pa 


B_  dx 
m 


B44  = P a0  U f b2  Bm  dx+p  ^CD b2  Bm  dx 


*u/dtfdx+,±CD /*,„ 


d d~f  dx 


B4*2=PaoU  J'dd2 -P3^CD/ y,0dd2  dx 

B46  =PaoU  J~xdd2dx+  P TrC»  f ^ 1 o x d d2  dx 
^46  = " A46/cj2  = "Pao  1)2  y dd2dx 

B5  5 =PaoU yx2Bm  dx+^CD y^lox2Bmdx 


TABLE  2 (Continued) 


B53  = "paoU /x  Bm  dx'p^CD  J dx 


C55 “ -pa 


„u2/«»„ 


dx 


0 u y*  x2ddx+p~cD  y*  y1(,x2ddx 


B66  P a 


B62  = p ao  U / x d dx  + p 


x d dx 


Fy*  = -p  oj, 


/xddx  +p  ^CD /vio 

B64=paoU  ylo  X dd2 

C66  =‘A66/w2  =-paoL'2y xddx 

f K (-d-,  + j x cos  (3) 

-0Aa0U y de  0 2 


dx 


cos  (K0  b(x ) sin  (3)dx 


ACd  / de 


/' 


Ko(-d2  +jxcos(3)  jK()bsin/3  ^ ( + e~j  K0 b sin  0 


ly,  J + c 


ly,pl>dx 


F->*  = -jp  to, 


/K  (-  d , + j x cos  (3 ) „ , 

Bm  e ° 1 cos  (K0  b sin  (3)  dx 


-JP 


4 r K (-  d.  + j x cos  (3 ) j K b sin  (3  -j  K b sin  (3 

^„*c BJ  Bmc  ' (e1  ° U|SI  + * I'.pl'dx 


r K (- d.  + j x cos  (3) 

F4*  = -pco0  Aa0U  J Bm  be  ° 1 sin  (KQ  b sin  (3)  dx 


4 r K (-d.  + j x cos  (3)  j K b sin  (3 

(e  u 


*1  s I - e 


-j  Ko  b sin  (3 


Iilp|)dx 


Fc*  = jp  oj 


r K J 

oAa0U J xBmc  ° 


- d . + j x cos  |3  ) 


cos  ( K()  b sin  (3)  dx 
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FIN- GENERATED  LIFT  COEFFICIENTS 

As  mentioned  earlier,  a SWATH  configuration  can  become  unstable  in  the  vertical-plane 
modes  at  high  speeds  due  to  an  inherently  small  waterplane  area.  Both  theoretical26  and 
experimental27  investigations  showed  that  small  horizontal  fins  fixed  on  the  inboard  sides  of 
the  hull  at  certain  longitudinal  positions  significantly  improve  the  stability  and  considerably 
reduce  the  peak  motion  amplitude. 

The  reduction  of  the  peak-motion  amplitude  foT  heave,  pitch,  and  roll  modes  was  more 
significant  when  the  stabilizing  fins  were  placed  both  forward  and  aft  of  the  longitudinal 
center  of  gravity.  For  stability,  the  aft  fins  should  be  larger  than  the  forward  fins.  The  mam 
factor  contributing  to  the  reduction  of  the  peak-motion  amplitude  stems  from  the  damping 
effects  of  the  fins.  It  will  be  shown  later  that  the  angle  of  attack  on  the  fins  is  proportional 
to  the  vertical  velocity  of  the  fins  when  the  ship  has  a forward  speed. 

Fins,  being  a lifting  surface,  generate  a lift  when  subjected  to  an  angle  of  attack.  For 
fins  attached  to  a SWATH  ship,  the  angle  of  attack  on  the  fins  can  be  generated  by  the  trim 
and  the  relative  fluid  velocity  with  respect  to  the  fins  in  the  horizontal  and  vertical  directions. 
The  trim  angle  can  easily  be  found  if  the  pitch  angle  of  the  body  is  known,  but  to  find  the 
relative  fluid  velocity  components  involves  the  motion  of  the  fins  with  respect  to  the  fluid 
motion  caused  by  the  incoming  and  diffracted  waves  and  the  forward  speed  of  the  ship.  It  is 
therefore  a tremendously  difficult  task  to  perform  an  accurate  evaluation  of  the  lift  generated 
by  the  fins.  To  name  a few  hydrodynamic  effects  to  be  considered,  there  are  body-fin  inter- 

“6L.ee.  C.M.  and  M.  Martin,  “Determination  of  Size  of  Stabilizing  Fins  for  Small  Waterplane  Area.  Twin-Hull 
Ships.”  DTNSRDC  Report  4495  (1974). 

'7Kallio,  J.A.  and  J.J.  Ricci,  “Seaworthiness  Characteristics  of  a Small  Waterplane  Area  Twin-Hull 
(SWATH  IV)  Part  II,”  DTNSRDC  Report  SPD  620-02  (1976). 
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actions,  the  blockage  effect  of  the  other  hull,  unsteady  effect,  boundary-layer  effect,  down- 
wash  and  upwash  effects  between  the  fore-and-aft  fins  and  from  the  fins  on  the  other  hull, 
the  f re  e-surface  effect  on  the  fins,  and  so  forth.  It  is  an  almost  impossible  task  to  take  into 
account  the  aforementioned  effects  satisfactorily,  hence  the  following  assumptions  will  be 
made  to  keep  the  analysis  in  a tractable  form.  Except  for  the  body-fin  effect  based  on  the 
investigation  performed  by  Pitts,  Nielsen,  and  Kaattari,28  all  other  effects  will  be  neglected. 
Justification  for  neglecting  some  of  the  effects  may  be  made,  but  most  of  the  justification 
will  be  based  on  heuristic  physical  arguments,  and  will  be  omitted  in  this  report  except  for 
the  investigation  shown  in  Appendix  B. 

Let  the  center  of  pressure  of  a thin  foil  of  span  S and  chord  C be  located  at  x = C.  The 
lift  L on  the  foil  with  a small  angle  of  attack  can  be  expressed  by 


L(t)  = |u2  A(f>CLaa(t)  (38) 

where  A*'  * is  the  plane  area  of  the  foil  and  CLQ,  the  lift-curve  slope.  The  C'La  will  be  assumed 
to  be  the  sum  of  the  two  effects  which  can  be  expressed  as 


where 


( LCK~  + ^BtWt'^LCX* 


w 


K 


iCLQ,) 

w,B,=lw 


W(B) 


W 


(39) 


K 


BfW) 


(Ct-a)BOV) 

^La^w 


in  which  the  subscript  W represents  the  case  of  the  foil  alone,  B(W)  the  case  of  the  lift  on 
the  body  induced  by  the  foil,  and  W(B)  the  case  of  the  lift  on  the  wing  induced  by  the  body. 
According  to  slender  body  theory,29,30  KW(B)  and  KB(W)  can  be  expressed  in  terms  of  the 
ratio  of  the  radius  of  the  body  r at  which  the  fin  is  attached  and  the  transverse  distance  from 
the  body  axis  to  the  tip  of  the  fin  rQ,  i.e., 


28Pitts,  W.C.,  J.N.  Nielsen  and  G.E.  Kaattari,  “Lift  and  Center  of  Pressure  of  Wing- Body-Tail  Combinations  at 
Subsonic,  Transonic,  and  Supersonic  Speeds,”  NACA  Report  1307  (1959). 

~9Spreiter.  J.R.,  “The  Aerodynamic  Forces  on  Slender  Plane-  and  Cruciform-Wing  and  Body  Combinations,” 
NACA  Report  962  (1950). 

'"Morikawa  G„  “Supersonic  Wing-Body  Lift,”  J.  Aeron.  Sci.,  Vol.  18,  No.  4,  pp.  217-228  (1951). 
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(40) 


(4!) 


where  6 = r/r().  The  foregoing  two  functions  are  plotted  in  Figure  6. 

An  empirical  formula  derived  by  Whicker  and  Fehlner31  for  low  aspect  ratio  wings  will 
be  used  for  the  evaluation  of  (CLa)w  The  expression  for  the  wings  of  zero  sweep  angle  is 
given  by 


1.8  rr  Ae 

(C’La)  = .....  per  radian 

1 .8  + •</  Ae“  + 4 


(42) 


Rather  than  taking  twice  the  value  of  the  geometric  aspect  ratio  as  did  Whicker  and  Fehlner. 
we  choose  the  effective  aspect  ratio  in  this  report  as 

j/C  (43) 


where  C is  the  average  chord  of  the  fin. 

The  numerator  of  the  above  equation  corresponds  to  the  distance  between  a point 
vortex  located  at  the  tip  of  the  wing  and  its  reflected  image  inside  a circle.  The  foregoing 
approach  is  mainly  designed  to  make  an  engineering  approximation  rather  than  seeking  a 
rigorous  rational  approach.  For  more  rigorous  approaches,  see  the  extensive  review  given  in 
Ashley  and  Rodden.32  Since  the  fins  are  attached  to  an  oscillating  body,  the  unsteady  effect 
on  the  lifting  characteristics  may  be  important.  However,  it  can  be  found  from  Lawrence 
and  Gerber33  that  the  unsteady  effect  on  (CLa)w  is  negligible  for  wings  of  aspect  ratio  less 
than  1.5. 

''Whicker.  L.F.  and  L.F.  Fehlner,  “Free-Stream  Characteristics  of  a Family  of  Low-Aspect-Ratio.  All- 
Movable  Control  Surfaces  for  Application  to  Ship  Design,”  David  Taylor  Model  Basin  Report  933  (1958). 

32  Ashley.  H.  and  W.P.  Rodden,  “Wing-Body  Aerodynamic  Interaction.”  Ann.  Rev.  Fluid  Mech..  Vol.  4 (1972). 

"Lawrence,  II. R.  and  F..H.  Gerber,  “The  Aerodynamic  Forces  on  Low  Aspect  Ratio  Wings  Oscillating  in  an 
Incompressible  Flow,"  J.  Aeron.  Sci.,  Vol.  19,  No.  1 1,  pp.  769-781  (1952). 
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The  angle  of  attack  on  a fin  at  the  center  of  pressure  located  at  (6,  bj , — d j CC ))  can  be 
obtained  under  the  assumption  of  a small  angle  of  attack  by 

(i3(t)-e£5(t)  + b,  i4(t)-fv(S,b1,-d.,t)) 
a(t)=£s(t)+  1 (44) 

where  the  wave-induced  vertical  velocity  f is  defined  in  Equation  (35).  The  center  of  pressure 
of  a fin  will  be  assumed  to  be  located  at  the  midspan  and  quarter-chord  from  the  leading  edge. 
Here  again,  the  wave-diffraction  effects  by  the  body  and  the  fins  are  neglected. 

When  a foil  oscillates  vertically,  a cross-flow  drag  and  virtual  inertia  force  result.  The 
cross-flow  drag  D(t)  can  be  expressed  as 


D(t)  = ^ A^f>CD(f>  (S3  - C£5  + bj  «4-fvm3-e£5  +b,  j4-fvl  (45) 

where  is  the  cross-flow  drag  coefficient  of  the  foil,  and  the  virtual  inertia  force  I(t)  can 
be  expressed  as 

I(t ) = ( m(f)  + a33<f>)  ($3  - C f 5 + b,  ) (46) 

Unless  the  mass  of  the  fin  is  specified,  it  will  be  approximated  by  taking  the  fin  section 
as  a neutrally  buoyant  ellipse  as  m*f*  = p 7r/4  SCt  where  S is  the  span  and  t is  the  maximum 
thickness  of  the  fin.  The  added  mass  will  be  approximated  by  a-j-^  = pS  n/4  C2.  For  N 
number  of  fins  for  each  hull,  the  heave  force  Fv*’\  pitch  moment  and  roll  moment 

contributed  by  2N  fins  can  be  obtained  from  the  foregoing  results  as  shown  below. 


Fv(f>  = l(t)  + L(t)  + D(t) 


= Z[(mi(f)  + a33i(f))(«3-8i^5> 
i = 1 L 

♦I  V2  A/f)  Ct  J 2 U * 2 ^ } 


Ai<0 4?{«3  - !i l's -fv<bli»  I b - «l i5  - i>l l> I 
+ «3-8i*5-M-bli))l53-*i^5-M-bli>l} 


(47) 


where  the  abbreviation  fy(b,j)  stands  for  fv(Cj.  bf  j,  - d,  j ) and  the  lift-curve  slope  for  the  ith 
fin  C,  aj  is  obtained  from  Equation  (39). 
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It  can  he  noticed  in  the  above  equation  that  the  roll-contributed  lift  and  drag  do  not 
appear  because  these  forces  on  the  fins  are  located  opposite  to  one  another  and  will  be 
canceled.  Strictly  speaking,  however,  the  vertical  velocity  contributed  by  the  roll  motion 
should  have  been  included  within  the  absolute  signs  in  Equation  (47).  Since  the  introduction 
of  the  roll  motion  in  the  vertical-plane  modes  complicates  the  computational  procedure,  it 
will  be  tacitly  assumed  that  the  roll  effect  on  Fv<*^  is  negligible.  The  pitch  moment  is 
obtained  by 

N 

MP<0  = -£eifi  <48 ) 

i = I 


where  fj  is  the  expression  given  within  the  brackets  after  the  summation  sign  in  Equation  (47). 
The  roll  moment  is  obtained  by 


MR(f)  = 


(m 


(f) 


+ a 


(f) 


3 3 i 


>$4 


+ ^ U Ai(f)CLai  {bIi(2b|jf4  -rv(bn)  + ?v(-bn))J 

+ 3Ai(f)CDi(f){blj(b1^4-fv(b|i».|bli^4-fv(bll,| 
+ bli(bii?4+fv(-bIi))|bljj4+fv(-bIi)|J 


(4‘M 


Table  3 indicates  hydrodynamic  coefficients  associated  with  the  stabilizing  fins.  The 
superscript  (I)  is  used  to  denote  these  coefficients.  Equilinearization  for  viscous  damping,  as 
shown  in  Equation  (37),  is  also  made.  The  relative  fluid  vertical  velocity  amplitude  at  the 
i th  fin  on  the  starboard  side  hull  zSQ  and  on  the  port  side  hull  z are  defined  by 


Zso=l$3-M5-M8i--bli--dli>l 

V = l$3-M5-M«i  bli'-dli>l  <51) 

and  the  relative  velocity  amplitude  due  to  the  roll  motion  is  also  similarly  defined  by 

ZsoR,  = lbti*4  + fv<8i--bli--dli>l  <52) 

Zpo(R,  = lb1i*4-fv<6i-bij.-dij)l  (53) 


38 


TABLE  3 - HYDRODYNAMIC  COEFFICIENTS  OF  STABILIZING  FINS 


i = 1 


A3/f>--H1(nfr>t,33p 

b,/' » p i A/f > { u cLoi  ^ cDim  ,iso  * ipo , | 

B3!'f ) - p I e,  A/')  | U Clal  * £ Cw<0  (zs0  + ipo  >} 
Cj5<n*pU2  I \lf,CLal 
A44°‘  1 bU2<m!ri*a33in> 

B44<n  * P I !>,  r A,*''  | « q.,  + £ CDff>  : » + ip.'* 1 

As5lfl*  1 *i2  (mff> + a33ff>> 

A53<fl " A35<r) 

b55<0  - p z t2  a p | u rLai  ♦ ± cDm  (zJ0  + ipo .) 
B53(f)-B35m 


C55<f).-pU2  I 6,A<f)CLoi 
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! 1/  L • O ■ >r  • rt 
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usn  c 


. -jK  b.-sinl3 
+ Zpoe  > 


F4(f>  = -pw0AUS  bn  A/0  CLaj  eK»('dli  + J C°S  P)  sin  (K0  b,  j si 
+ Jp^«oAZbI1A/0CDf0eV--,.*J«,«-« 

(i  (R)  ejKobliSin^_  • (R)  "j  Koblisin^v 

'so  'po  L ’ 


->0  A u Z b,  J AM)  CLai  e 0 “>*  J *•  ^ sin  (K0  b, , sin  0) 

-a;  A I b,.A.(0C  (O  eKo(-d,i  + j Ci cos^ 
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TABLE  3 — (Continued) 


FS<0  - j p u,  A U I «,  A/' , CL„f  ,K"'-d"  * 1 'l « cos  , Kn  b, , sin  (J, 


*jPi;co0Ari!iA(n(:D1(f)c'!.(-d.,*i«,^« 


<Z, 


. ejKob„sin0+  . c“i  KQ  b(j  sin  (J 


po 


The  hydrodynamic  coefficients  of  a SWATH  ship  with  stabilizing  fins  are  the  sum  of 
the  effects  of  the  potential  How,  viscosity,  and  fins;  the  final  expressions  for  the  hydrodyna- 
mic coefficients  to  be  used  in  the  equations  of  motion  ((5)  and  (6))  can  be  obtained  by 
summing  the  coefficients  in  Tables  I through  3 and  are  shown  in  Table  4.  An  exception  to 
the  straightforward  summation  is  made  for  B35,  BS3.  and  A5S  in  which  A33  is  taken  as  the 
sum  of  the  heave  added  mass  of  the  bare  hull  and  the  fins. 


TABLE  4 - TOTAL  HYDRODYNAMIC  COEFFICIENTS 
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TABLE  4 - (Continued) 
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TABLE  4 (Continued) 
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TABLE  4 (Continued) 
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TABLE  4 (Continued) 
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TABLE  4 (Continued) 
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RESULTS  AND  DISCUSSION 


Model  experimental  results  available  at  DTNSRDC  are  now  utilized  to  venfy  the  validiiy 
of  the  theoretically  predicted  values.  The  comparisons  are  made  for  two-dimensional  hydro- 
dynamic  coefficients  of  twin  cylinders  of  semisubmersible  cross  section  as  well  as  for  three- 
dimensional  hydrodynamic  coefficients  of  the  single  hull  of  a SWATH  model  and  of  complete 
twin  hulls  with  and  without  fins.  The  motion  of  two  SWATH  models,  SWATH  4 and  SWATH 
OA,  in  angular  waves  are  also  compared. 

The  particular  dimensions  of  the  SWATH  models  used  to  obtain  experimental  data  are  . 

1 

shown  in  Table  5 along  with  the  sizes  and  locations  ot  the  stabilizing  fins  attached  to  SWATH 
6 A. 

More  complete  descriptions  of  the  experiments  for  the  models  shown  in  Table  5 have 
already  been  published.^7,14'* 

Since  the  main  backbone  of  the  strip  theory  is  based  on  two-dimensional  hydrodynamic 
coefficients,  it  is  important  to  check  the  validity  of  the  two-dimensional  hydrodynamic 
calculations. 

At  an  early  stage  in  the  SWATH  program,  a forced  oscillation  of  twin  cylindrical  bodies 
of  semisubmerged  cross  section  was  conducted  at  the  Center  by  D M.  Gerzina.  The  forced 
oscillation  was  made  independently  in  heave,  sway,  and  roll  mode.  The  dimensions  of  the 
cross  section  of  the  model  and  the  experimental  results  are  presented  in  Figures  7-9.  The 
theoretical  curves  presented  in  these  figures  were  obtained  by  the  method  of  source  distribu- 
tion.71 Two  or  three  amplitudes  of  oscillation  were  used  in  the  experiment  to  check  the 
linearity  of  the  results.  The  amplitudes  shown  in  Figure  9 are  the  vertical  strokes  of  the 
pivoted  arms  used  to  generate  the  roll  oscillation.  The  notation  a (=  10.15  cm)  is  the  radius 
of  one  hull,  to  is  the  circular  frequency  of  oscillation.  SA  is  the  submerged  sectional  area  of 
the  two  hulls,  and  b (=  30.48  cm)  is  one-half  the  distance  between  the  centerline  of  each  hull. 

The  distinct  discontinuities  apparent  at  certain  frequencies  in  the  figures  are  due  to  the 
trapped  mode  of  standing  waves  between  the  two  bodies.  Although  it  appears  that  heave 
oscillation  did  not  include  discontinuities  in  the  hydrodynamic  coefficients,  it  would  occur  at 
the  higher  frequencies.  The  sway  and  heave  added  mass  coefficients  (a-,-,  and  a A A ) obtained 
theoretically  and  experimentally  correlate  fairly  well.  The  discrepancy  observed  in  Figure  9 


14Lee,  C M.  and  L.O.  Murray,  “Experimental  Investigation  of  Hydrodynamic  Coefficients  of  a Small  Water- 
plane  Area.  Twin-Hull  Model  (SWATH  6AI.”  DTNSRDC  Report  (in  preparation). 

•Reported  informally  for  SWATH  I by  R.  Stahl. 
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TABLE  5 - PARTICULARS  OF  THE  SWATH  MODELS 
USED  TO  OBTAIN  DTNSRDC  EXPERIMENTAL  DATA 

(Dimensions  are  full  scale) 


Parameter 

SWATH  1 

SWATH  4 

SWATH  6A 

Displacement,  long  tons  in  s.w. 

2S543 

41  15 

2802 

Length  at  the  Waterline,  m 

1 34 

69. 1 

52.5 

Length  of  Main  Hull,  m 

158.5 

87.7 

73.2 

Beam  of  Each  Hull  at  the  Waterline,  m 

5.2 

2.4 

Hull  Spacing  between  the  Centerlines,  m 

42. 1 

22  9 

22.9 

Draft  at  the  Midship,  m 

17.3 

9.1 

8.1 

Bridging  Structure  Clearance  from 
Waterline,  m 

5.2 

6.1 

Maximum  Diameter  of  Main  Hull,  m 

9.4 

5.5 

4.6 

Longitudinal  Center  of  Gravity  Aft  of 
Main  Hull  Nose,  m 

77.5 

42.3 

35.5 

Vertical  Center  of  Gravity  above 
Baseline,  m 

15.3 

9.6 

10.4 

Transverse  GM.  m 

- 

2.4 

2.9 

l ongitudinal  GM.  m 

- 

1 1.1 

6.8 

Radius  of  Gyration  For  Pitch 

3 2.8 

20.8 

16.9 

Radius  of  Gyration  For  Roll 

- 

- 

10.2 

Waterplane  Area,  m 

1272.7 

251.3 

193.9 

Model  Scale  Ratio 

1/40.94 

1/20.4 

''”5 

SWATH  6A  Fins  (Rectangular  Planform): 


Chord 

Span 

Location 

m 

m 

m 

Fore  Fin 

2.59 

3.1  1 

1 7. 1 5 

Aft  Fin 

4.48 

5.36 

62. 24 

•Distance  from  the  main  hull  nose  to  the  quarter-chord  point. 

for  the  roll  added  inertia  a44  is  construed  as  resulting  from  experimental  errors  due  to 
difficulty  in  generating  pure  roll  oscillation.  The  sway  and  roll  damping  coefficients  (b-,-,  and 
b44 ) showed  considerable  discrepancies  between  theory  and  experiment.  The  theoretical  values 
of  damping  coefficients  were  not  expected  to  correlate  well  with  the  experimental  results 
because  of  the  limitation  of  the  potential-flow  assumption,  especially  for  the  roll  damping 
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Figure  8 - Heave  Added  Mass  and  Damping  Coefficients  of 
Semisubmersible  Twin  Cylinders 
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Figure  9 - Roll  Added  Inertia  and  Damping  Coefficients  of 
Semisubmersible  Twin  Cylinders 
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coefficient.  The  heave  damping  coefficient  (b^)  in  Figure  8 showed  good  agreement  between 
theory  and  experiment,  except  for  higher  frequencies.  Similar  experimental  results  for  b,-,  in 
a single  cylinder  of  bulbous  cross  section  in  low  frequencies  showed  a distinct  discrepancy  with 
theoretical  values,  as  shown  in  Figure  5 in  which  the  dividing  factor  for  b3 3 
was  the  same  as  SA.  The  twin  cylinders  can  he  expected  to  exhibit  a discrepancy  similar  to 
that  observed  for  the  single  cylinder  in  the  small  frequency  region.  In  fact,  this  will  be 
demonstrated  later  when  the  three-dimensional  heave  damping  coefficient  is  discussed. 
Unfortunately,  the  frequency  range  chosen  in  the  experiment  for  the  twin  semisuhmersible 
cylinders  was  rather  high  from  the  viewpoint  of  SWATH  motion  The  horizontal  scales  in 
Figure  7-9  are  equivalent  to 

2 rra/X  (=  a g ) 

where  X is  the  wave  length. 

For  oj~  a/g  = 0.2, 

X = 1 0 7r  a 

which  is  about  the  magnitude  of  one  ship  length  of  a SWATH  ship.  As  can  be  seen  later  in 
the  comparisons  of  motions,  the  SWATH  ships  hitherto  examined  were  barely  excited  by  a 
wave  of  such  short  length. 

Comparisons  of  three-dimensional  hydrodynamic  coefficients  are  shown  in  Figua-s  10 
through  12.  Figure  10  compares  the  theoretical  and  experimental  values  of  heave  added  mass 
coefficient  A33  and  damping  coefficient  of  SWATH  1.  Both  demihull  and  twin-hull 
results  are  presented  for  A33.  It  is  interesting  to  note  that  the  demihull  theory  appears  to  fit 
the  twin-hull  experimental  results  well.  The  dip  exhibited  by  the  twin-hull  theory  may  appear 
suspicious:  however,  such  a dip  in  A33  has  been  confirmed  for  a conventional  catamaran  by 
another  experiment.6  The  dip  in  A^  for  a twin-hull  ship  is  considered  to  be  caused  by  a 
mutual  blockage  effect  between  the  two  hulls.  Comparison  of  the  theoretically  predicted  dip 
in  A33  for  a SWATH  shape  cannot  be  made  at  present  due  to  the  lack  of  available  experimental 
results.  The  theoretical  values  of  A33  for  co^/T/g  < 2.0  appeared  to  be  lower  than  the 
experimental  values. 

As  demonstrated  earlier  in  Figures  3 and  4,  the  damping  coefficients  involved  in  the 
vertical-plane  modes  play  a sensitive  role  in  the  prediction  of  heave  and  pitch  motion.  As 
pointed  out.  the  damping  coefficients  obtained  under  potential-now  theory  for  small-amplitude 
oscillation  yield  overestimated  peak  amplitudes  of  the  motion  at  the  resonant  frequency 
compared  to  the  magnitude  measured  in  motion  experiments  in  regular  waves.  The  wave 
amplitudes  used  in  the  motion  experiments  were  such  that  the  ratio  of  the  wavelength  to  the 
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amplitude  was  usually  greater  than  100;  this  can  normally  be  coasidered  to  be  adequate  for 
application  in  the  development  of  a linear  boundary-value  problem  involving  a free  surface. 
Unfortunately,  this  turned  out  to  be  untrue  in  the  prediction  of  SWATH  motion  at  the 
resonant  frequencies. 

Figure  10  presents  three  theoretical  values  of  heave  damping  coefficients,  one  for  the 
demihull  and  the  other  two  for  the  twin  hulls.  One  of  these  two  was  obtained  under  the 
potential-flow  assumption  which  is  expressed  in  Table  I as  B33  and  the  other  included  the 
so-called  cross-flow  effects,  denoted  by  in  Table  2.  Theoretical  values  for  the  demihull 
seemed  to  agree  well  with  the  experimental  results  for  both  the  demihull  and  twin  hulls 
obtained  under  small  amplitudes  of  oscillation.  However,  the  demihull  theory  appeared  to 
underprediet  damping  in  the  small  frequency  region,  cj  yj  L/g  < 2.0,  for  large- amplitude 
oscillation.  Note  that  even  the  largest  amplitude  of  oscillation  (5.1  cm)  corresponds  to  only 
about  one-tenth  the  draft  of  the  model.  At  the  resonant  frequency  the  heave  amplitude  of  a 
SWATH  model  without  fins  can  reach  about  three  times  the  wave  amplitude,  or  about  15  cm 
for  the  SWATH  I model  at  zero  speed  if  the  ratio  of  the  wavelength  to  amplitude  is  assumed 
to  be  100.  Thus,  the  amplitude  of  5.1  cm  does  not  reflect  the  peak  amplitude  of  oscillation 
of  a SWATH  model  in  waves.  As  can  be  seen  in  Figure  10.  the  heave  damping  coefficients 
obtained  by  the  strip  theory  under  the  potential- flow  assumptions  were  considerably  lower 
than  the  experimental  values  in  the  high  frequency  range.  The  chained  curve  which  represents 
the  present  theory  including  the  cross-flow  effects  seems  to  fit  the  experimental  resuits  better 
at  high  frequencies. 

It  is  of  more  interest  to  us  to  examine  the  damping  coefficients  of  a SWATH  ship  when 
it  has  a forward  speed.  Furthermore,  as  mentioned  earlier,  SWATH  ships  should  have 
stabilizing  fins  for  stability  in  high-speed  operations.  The  stabilizing  fins  can  contribute 
significant  damping  in  heave,  pitch,  and  roll  modes  in  high-speed  range.  Figure  1 1 compares 
the  heave  and  pitch  damping  coefficients  obtained  by  the  present  theory  and  experiments. 
Three  hull  conditions  were  examined:  the  bare  hull,  the  hull  with  aft  fins,  and  the  hull  with 
fins  fore  and  aft.  Table  5 indicates  sizes  and  locations  of  the  fins  used  in  the  experiment. 

The  fins  were  attached  on  the  inboard  side  of  each  hull  and  were  stationary.  The  theoretical 
results  shown  were  obtained  by  the  expressions  given  in  Table  4.  The  values  of  coefficients 
chosen  were  as  follows:  aQ  = 0.07,  CD  = 0.5,  C'La  for  the  fore  fin  = 4.38.  and  Cla  for  the 
aft  fin  = 3.2.  The  C(  a values  were  obtained  by  using  Equations  (39).  (42).  (43)  and  Figure  6. 
The  theoretical  and  experimental  values  agreed  very  well.  The  bar  signs  designate  nondimen- 
sional  quantities  which  are  define  as 
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where  L is  the  main  hull  length.  Comparisons  of  other  coefficients  associated  with  the  heave- 
pitch  motion  are  given  in  Lee  and  Murray.34  Most  of  the  trends  were  in  good  agreement  for 
the  frequency  of  oscillation,  hut  some  of  them  indicated  almost  constant  magnitudes  of 
difference  between  the  two  results  over  the  range  of  frequencies  examined  in  the  experiment. 
Since  there  have  been  some  unresolved  uncertainties  in  the  experimental  measurements,  it  is 
difficult  to  judge  at  present  whether  theoretical  prediction  of  those  coefficients  is  unreliable. 

At  any  rate,  uncoupled  damping  coefficients  constitute  the  most  sensitive  coefficients  with 
respect  to  peak  amplitudes.  The  good  agreement  obtained  for  B33  and  B35  is  the  major 
con tnbu ting  factor  that  produced  good  agreement  in  the  vertical  plane  motion  shown  in  the 
subsequent  figures. 

Figure  12  shows  the  wave-exciting  heave  amplitudes  and  phases  of  SWATH  6A  with  fore- 
and-aft  fins  versus  the  wave-toship  length  ratio.  The  ship  length  is  that  of  the  main  hull.  i.e.. 
73.15  m.  Agreement  between  theoretical  and  experimental  results  was  good  except  the  phase 
angles  at  zero  speed.  Similar  comparisons  for  wave-exciting  pitch  moment  were  poor.  How- 
ever. for  the  same  reasons  stated  previously,  it  is  not  obvious  whether  the  theoretical 
prediction  is  poor.* 

Comparisons  of  motion  are  shown  in  the  next  seven  figures.  Figure  13  shows  the 
amplitudes  and  phases  of  the  heave  and  pitch  motion  of  SWATH  4 without  fins  in  regular 
head  waves  at  20  knots.  The  results  are  given  versus  the  frequency  ratio  oj/cJn  where 
con(=  0.6128  rad/sec)  is  the  natural  heave  frequency.  It  can  be  seen  clearly  from  Figure  13 
that  SWATH  ships  exhibit  a highly  tuned  response  at  the  heave  natural  frequency.  This  highly 
tuned  motion  of  SWATH  ships  does  not  usually  occur  for  monohull  surface  ships,  except  for 
roll  motion.  The  results  presented  in  Figure  13  demonstrate  the  validity  of  the  theoretical 
prediction  of  motion  for  the  bare-hull  condition.  The  effects  of  the  viscous  coefficients  aQ 
and  CD  in  reducing  the  peak  amplitudes  can  be  seen  by  comparing  Figures  3 and  13. 

Figure  14a  presents  the  amplitudes  of  heave,  pitch,  and  relative  vertical  motion  at  the 
edge  of  the  forward  deck  of  SWATH  6A  with  fore-and-aft  fins  in  regular  head  waves  at  0 and 
28  knots.  An  accurate  position  of  the  relative  motion  was  2.74  m from  the  main  hull  nose. 


i 


a 


•The  experimental  values  were  consistently  almost  twice  the  theoretical  values  of  amplitude  over  the  range  of 
wavelengths.  If  the  theoretical  values  were  wrong,  the  agreements  shown  in  Figures  13,  14a,  and  14b  for 
pitch  amplitudes  and  phases  could  not  have  been  obtained. 
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Normally,  the  relative  motion  prediction  is  less  reliable  than  other  motions  because  the  theory 
does  not  account  lor  the  deformation  of  the  incoming  waves  caused  by  the  diffraction  by  the 
body  and  by  the  body  motion.  However,  agreement  between  the  two  results  was  very  good 
in  this  case,  particularly  at  28  knots.  It  suggests  that  the  wave  deformation  near  the  bow  may 
not  be  so  significant  as  in  the  case  of  conventional  ships. 

Motion  amplitudes  of  SWATH  6A  at  20  knots  are  respectively  presented  in  Figures  14b 
and  14c  for  bow-quartering  waves  and  following  waves.  Agreement  between  the  two  results  is 
good  in  the  case  of  bow-quartering  waves.  The  largest  discrepancy  in  pitch  amplitude  between 
the  wave-to-ship  length  ratios  of  3 and  4 was  in  the  order  of  experimental  error,  i.e.,  about  ±0.03  deg. 

As  expected,  the  prediction  of  motion  in  following  waves  was  poor.  At  20  knots,  the 
wavelength  which  gave  zero  encountering  frequency  was  0.93  times  the  length  of  the  ship. 

When  the  encountering  frequency  is  extremely  small,  the  hydrodynamic  coefficients  obtained 
by  the  strip  theory  becomes  invalid  since  the  strip  theory  is  based  on  the  assumption  of  high 
frequency  of  oscillation.20  The  large  amplitudes  predicted  by  the  theory  near  X/L=  1.0  are 
considered  to  be  errors  resulting  from  the  shortcomings  of  the  strip  theory  for  small  encoun- 
tering frequencies.  Furthermore,  as  observed  in  the  experiment,  a SWATH  model  exhibits 
some  instabilities  in  the  vertical-plane  mode  near  the  zero  encountering  frequency.  Thus  it 
can  easily  be  judged  that  any  theory  based  on  harmonic  time  dependence  would  fail  to  pro- 
vide a realistic  solution.  The  experimental  results  shown  in  Figure  14c  should  by  no  means 
be  interpreted  to  be  as  reliable  as  those  for  head  waves.  For  those  following  waves  which  can 
make  the  encountering  frequency  small,  it  is  extremely  difficult  to  keep  a model  at  constant 
speed:  hence  an  effort  by  the  person  who  controls  the  propeller  rpm  to  maintain  a constant 
speed  of  the  model  can  impose  a significant  man-made  surge  motion  on  the  model.  The  large 
surge  motion  generated  in  such  a manner  seems  to  induce  motions  in  other  modes. 

Although  the  present  method  of  predicting  the  motion  of  a SWATH  ship  in  following 
and  stem-quartering  waves  may  still  be  useful,  it  definitely  needs  some  improvement. 

Figure  15  shows  the  roll  amplitude  of  SWATH  6A  at  0 and  20  knots  in  regular  beam 
waves.  Agreement  between  theoretical  and  experimental  results  is  good  except  for  the  peak 
amplitude  at  the  roll  resonance  at  zero  speed.  It  appears  that  larger  values  of  CD  and  CDi<f> 
than  those  used  for  heave  and  pitch  motion  may  be  necessary  to  reduce  the  peak  amplimde. 

Except  for  the  roll  resonant  frequency  at  zero  speed,  roll  amplitudes  were  less  than  1 deg  for 
A A > 100.  The  highly  tuned  roll  response  at  the  zero  speed  disappeared  at  20  knots  and  this 
was  correctly  predicted  by  the  theory.  The  damping  contributed  by  the  fins  and  the  viscous 
effects  seem  to  be  the  main  factors  which  reduced  the  resonant  roll  amplitudes. 


58 


figure  14.i 
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Figure  16  indicates  die  correlation  of  motion  of  a geosim  of  SWATM  6A  with  a displace- 
ment of  14.100  long  tons  in  irregular  waves.  The  irregular  waves  used  for  the  computation 
were  assumed  to  be  unidirectional  and  propagating  opposite  to  the  ship  heading.  The  wave 
spectra  obtained  at  Station  India  in  the  North  Atlantic  Ocean13  were  used  with  the  proper 
weighting  factor  for  their  frequency  of  occurrence. 14  The  RAO’s  obtained  by  the  present 
theory  were  used  together  with  the  aforementioned  spectra  to  obtain  the  relative  bow  motion 
represented  by  the  square  symbols  in  the  figure.  The  equation  to  obtain  these  values  is  given 
by 


( R) 


S(wo)doJo 


1/2 


(54) 


where  S(cj0>  is  tire  wave  spectrum,  and  £V,R)  = £V|R)/A  is  obtained  by  Equation  (II).  The 
solid  line  in  the  figure  was  obtained  by  using  the  Pierson-Moskowitz  formula  (see  liquation 
(16)  for  S(ce0)i  in  Equation  (54),  and  the  circle  mark  was  obtained  by  using  the  wave  spectrum 
employed  for  the  towing  basin  experiment.  The  triangular  mark  was  obtained  entirely  from 
the  experiment.  The  agreement  between  experimental  and  theoretical  values  when  an  identical 
spectrum  is  used  is  very  good.  Data  obtained  by  the  Station  India  spectra  and  the 
Pierson-Moskowitz  spectrum  formula  are  included  merely  to  indicate  that  motion  data  obtained 
by  using  the  basin-generated  irregular  waves  are  close  to  those  obtained  by  the  other  two 
methods. 

Figure  17  shows  (he  probability  that  the  significant  values  and  the  most  probable 
extreme  values  of  the  heave  amplitude  of  the  14.100-ton  SWATH  6A  will  be  exceeded  during 
operation  in  the  North  Atlantic  Ocean.  The  most  probable  extreme  values  were  obtained 
from  Equation  (19)  by  letting  6=  I and  T=68-  1.4  Hs  where  H$  is  the  significant  waveheight 
in  feet.  This  way  of  presenting  such  motion  results  as  slamming,  deck  wetness,  and  vertical 
accelerations  is  deemed  useful  in  a comparative  design  study  for  seakeeping  qualities  of  two  or 
more  candidate  hull  forms. 
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Figure  If)  Significant  Relative  Bow  Motion  Amplitudes  for  SWA  I II  (>A  in  Irregular 


PROBABILITY  OF  EXCEEDANCE  (PERCENT! 


Figure  1 7 Probability  of  SWATH  6A  Exceeding  the  Significant  and  Most  Probable 
Extreme  Values  of  Heave  Amplitude  at  26.3  Knots  in  the  North  Atlantic  Ocean 
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SUMMARY  AND  CONCLUSIONS 


This  report  presents  a theoretical  method  tor  predicting  the  motion  of  small-waterplane- 
area.  twin-hull  (SWATH)  ships  in  waves.  The  method  divides  equations  of  motion  into  three 
independent  groups;  equations  for  surge  motion  alone,  equations  for  coupled  heave  and  pitch 
motion,  and  equations  for  coupled  sway,  roll,  and  yaw  motion  Ihe  equations  are  based  on 
the  linear  frequency  response  of  a ship  to  harmonic  wave  excitation.  The  reference  frame  of 
the  measure  of  the  oscillatory  motion  of  a ship  is  a ( artesian  coordinate  system  translating 
on  the  calm-water  plane  parallel  to  the  mean  heading  of  the  ship  with  a constant  mean  speed. 
When  the  ship  is  at  rest  in  calm  water,  the  origin  of  the  coordinate  system  is  located  on  the 
free  surface  directly  above  or  below  the  ship  longitudinal  center  of  gravity. 

The  hydrodynamic  coefficients  associated  with  the  equations  of  motion  are  assumed  to 
consist  of  three  parts.  The  first  part  is  obtained  under  the  potential-flow'  assumption  (see 
Appendix  A).  The  basic  ingredient  of  this  part  comes  from  the  solution  of  two-dimensional 
boundary-value  problem  for  oscillating  twin  cylinders  in  a free  surface.17  The  expression  tor 
the  hydrodynamic  coefficients  contributed  by  the  first  part  is  given  in  Table  I. 

The  second  part  is  obtained  by  an  empirical  method  derived  from  the  cross-flow  approach 
to  a slender  body  at  a moderate  angle  of  incidence  in  a uniform  flow  74  A psuedo-steady- 
statc  assumption  is  invoked  to  apply  the  method  loan  oscillatory  body.  Two  coefficients 
involved  in  the  cross-flow  approach  are  the  viscous-lift  coefficient  a which  is  associated  with 
the  forward  velocity  of  the  ship  and  the  cross-flow  drag  coefficient  CD  which  is  independent 
of  the  forward  speed. 

An  equivalent  linearization  method  is  used  to  include  the  cross-flow'  drag  tenn  as  a linear 
damping  term  in  the  equations  of  motion.  This  necessitates  an  iterative  solution  of  the 
equations  of  motion  until  a reasonable  convergence  of  the  solution  is  obtained. 

The  third  part  is  that  contributed  by  the  stabilizing  fins,  if  any.  This  part  is  obtained  by 
using  the  method  developed  by  Pitts  et  al.’8  The  analysis  accounts  for  the  lift  induced  by  the 
fin  on  the  body  and  vice  versa. 

Theoretical  predictions  were  verified  by  comparing  them  with  available  results  from  model 
experiments.  The  sectional  added  mass  and  damping  coefficients  of  a SWATH  cross  section, 
the  three-dimensional  hydrodynamic  coefficients  associated  with  heave  and  pitch  motion,  and 
the  wave-exciting  heave  force  were  checked.  The  predicted  motions  of  two  SWATH  models 
were  compared  with  model  experimental  data  to  examine  the  effects  of  fins,  wave-heading 
angle,  ship  forward  velocity,  and  wavelength.  The  predicted  motion  in  irregular  waves  was 
also  compared  with  model  experimental  data. 


From  the  present  study,  the  following  conclusions  have  been  drawn 

1.  The  linear  frequency  response  theory  can  be  applied  for  predicting  motions  of 
SWATH  ships  in  waves. 

2.  The  strip  theory  used  in  deriving  the  hydrodynamic  coefficients  is  acceptable  a- 
an  approximation  in  predicting  motions  for  SWATH  ships. 

3.  The  effects  of  viscous  damping  on  SWATH  motions  can  be  predicted  by 
combining  the  cross-flow  approach  with  strip  theory.  In  a strict  theoretical  sense,  an 
extension  of  the  equilinearization  method  to  predict  motion  in  irregular  waves  (as  presently 
used)  cannot  be  justified;  however,  for  practical  purposes,  the  present  approach  is  recommended 
until  a more  rigorous  and  practical  means  is  developed. 

4.  The  present  method  of  predicting  the  effects  of  stabilizing  fins  on  the  motion  of 
a SWATH  ship  in  waves  seems  to  be  acceptable  despite  neglect  of  such  seemingly  important 
effects  as  free  surface,  unsteadiness,  and  downwash  on  the  aft  fins  by  the  forward  fins. 

5.  Prediction  of  SWATH  motion  in  following  waves  is  an  area  which  needs 
improvement  of  both  theory  and  model  experimental  techniques.  This  is  particularly 
important  if  future  plans  call  for  equipping  SWATH’s  with  controllable  fins  to  maintain  good 
seakeeping  quality  in  following  waves. 
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APPENDIX  A 

DETERMINATION  OF  HYDRODYNAMIC  COEFFICIENTS 
UNDER  POTENTIAL  FLOW  ASSUMPTION 


If  we  assume  that  the  fluid  surrounding  a catamaran  has  irrotational  motion,  we  can 
introduce  a velocity  potential  4mx,  y.  z.  t)  in  the  fluid  region.  The  velocity  potential  should 
satisfy,  in  addition  to  the  Laplace  equation,  the  following  boundary  conditions 


•iMx.  y.  z.  t)  + g = 0 on  z = 0 


where  g is  the  gravitational  acceleraton. 


V4>  • n = Vp  on  S0 


where  n is  the  unit  normal  vector  on  the  body  surface  pointing  into  the  body.  Vn  is  the 
normal  component  of  the  velocity  of  the  body  surface,  and  SD  is  the  mean  position  of  the 


<t>z(x.  y.  t)=  0 (57) 

^ i 1/2 

and  a physically  appropnate  far-field  condition  for  (x-  + y ~ ) ' -*'00  to  make  the  problem 
well  posed. 

Assuming  that  the  flow  disturbance  is  a small  perturbation  from  uniform  flow,  we  can 
express  in  the  form 

<!>  = - U x + 0s<x.  y.  z l + 0o<x.  y.  z ) e_l  1 58 ) 

Here  U is  the  forward  velocity  of  the  ship.  0 ^ is  the  steady  potential  representing  the  wave- 
making disturbance  due  to  the  forward  velocity,  and  0O  = 0OC  + j 0OS  is  the  potential 
associated  with  the  oscillatory  fluid  disturbance.  We  can  further  divide  0(1  with  three  distinct 
origins  of  the  oscillatory  fluid  disturbance  as 

0O  = 0,  + 0D  + 0M  <?l)> 

The  incoming  wave  potential  0,  is  given  by 

. _ jgA  K (z  + jx  cos0- jy  sin  0) 

0. e ((nil 

1 cu 

o 


Here  (3  is  the  wave-heading  angle  with  respect  to  the  positive  x-axis.  co()  is  the  wave  frequency. 
K(,  is  the  wave  number  given  by  co0"/g,  and  A is  the  wave  amplitude.  The  wave  diffraction  is 
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represented  by  0()  and  the  fluid  disturbance  caused  by  the  motion  of  the  body  in  initially 
calm  water  is  represented  by  0M  Within  a linear  approximation  to  the  solution  of  the  velocity 
potential,  we  can  let 


u 

^St  ^ 1 ^ko^k 


where  0k  is  another  set  of  velocity  potentials,  and  £ko,  k = 1,2 6 are  the  complex 

amplitudes  of  the  displacement  of  the  body  from  its  mean  position  in  surge,  sway,  heave,  roll, 
pitch,  and  yaw  modes,  respectively.  The  pressure  at  any  point  of  the  hull  is  obtained  from 
the  Bernoulli  equation  by 


^ U-  - p ^<l>t  + gz(t)  + I V4>  hj 


At  this  point,  we  will  establish  the  following  conditions.  ( 1 ) the  motion  of  the  body  is 
small  and  so  the  pressure  at  a point  on  the  body  surface  at  any  instant  can  be  obtained  via 
the  Taylor  expansion  of  the  pressure  at  the  mean  position  of  the  body:  (2)  the  terms 
Oi0s:.  0$0O.  0Q2.  £ko0s,  £ko0D.  £ko  ^ wil1  be  discarded  in  the  evaluation  of  pressure;  (3)  only  those 
components  of  the  pressure  which  have  harmonic  time  dependence  will  be  considered;  and 
(4)  the  static-pressure  component  and  the  component  which  contributes  to  the  static  restonng 
force  or  moment  of  the  body  will  not  be  included  in  the  evaluation  of  pressure. 

With  the  foregoing  conditions,  the  complex  amplitude  of  the  pressure  at  a point  on  the 
body  surface  can  be  expressed  by 


p = p jeo  + U 


= P (juJ  + U ~t)  (0|  + 0d  + *ko  ^k 


evaluated  at  the  mean  position  of  the  body. 

Integration  of  the  pressure  over  the  wetted  surface  of  a ship  should  yield  the  hydrody- 
namic forces  and  moments.  Thus,  we  have 


T"  = R 


pnjds 


JJn.  (jw  + u U + + *ko  M d" 


tor  i = I.  2.  3,  4.  5.  6.  where  n j , n2.  and  n3  are  respectively  the  x-,  y-,  and  z-components  of 
the  unit  normal  vector  and  n4  = yn3  - zn2,  n$  = znj  - xn3,  n&  = xn2  - yn^ 

We  can  turther  decompose  the  hydrodynamic  force  into  two  parts,  i.e., 

H>  _ p(e  ) . r-(m  ) 
i i i 

where 

p.(c)  = wave-excited  force 


and 


(0,  + 0D  ) ds 


(63) 


F.(m) 

= motion-excited  force 


-p  £ hoJJJ'^  + V ^)^kds  (64) 

k = 1 S0 

Applying  the  results  of  Ogilvie  and  Tuck20  to  the  simplified  case  in  this  work,  and  with 
the  use  of  the  Kronecker  delta  function  6^,  we  can  show  for  any  differentiable  scalar  function 
0 that 


dx 


n25i6 


) 0 ds 


-l 


C(x  ) 


ns  0 d£ 


(65) 


where  S is  the  immersed  hull  surface  forward  of  the  cross  section  at  x,  and  C(x)  is  the  line 
integral  along  the  contours  of  the  cross  section.  Utilization  of  Equation  (65)  in  Equation  (63) 
yields 


Fi(e,  = P JJ n{  (j^  + u|;)  (0,  +0D)  ds 
S0 

= P JJ  *n\  + U K0  cos  j3)0,  ds 

So 

+ P JJ + Un3  5jS  - U n2  6j6)0n  ds 


(66) 


where  the  line  integral  at  the  stem  section  is  neglected. 


Within  the  accuracy  of  the  order  of  approximation  in  this  work  and  together  with  a 
slenderness  assumption  tor  each  hull  of  the  ship,  i.e.,  rij  = o(n2  or  n3 ),  we  can  derive  from 
Equation  (56)  and  (59) 

30m 

— — = - joj  ( £ + a x r)  • n - U (a  x e , ) • n on  S (67) 

3n  _o  -o  - — -u  -i  - 

where 

Jo  - ^ 1 o'  ^2o’  ^3  o * 

-o  ~ ^4o’  ^5o' 
r=  (x,  y,  z) 

and  e,  is  the  unit  vector  in  the  x-direction.  The  original  derivation  of  Equation  (67)  without 
the  slenderness  assumption  was  derived  by  Tiinman  and  Newman.35  Substitution  of 
Equation  (61)  into  Equation  (67)  yields 

30j 

— = -jo;  n;  + U n3  6j5  - U n2  5i6  <68> 

on  SQ  for  i = 2.  3,  4,  5.  6.  From  Equation  (68),  we  can  deduce  the  following  relations: 


in  which  approximations  n5  =*  x n3  and  n6  — x n2  were  used. 
Using  Equation  (68)  and  the  relation 


co  + U KQ  cos  (3  = coc 


in  Equation  (66),  we  get 


F/e)  = jpojQ  ff  nj0,  ds 
J\ 

~ P ff  (*i"  + f ^3n  5iS  * f *2n  5i6)  * 


^Timman.  R.  and  J.N.  Newman,  “The  Coupled  Damping  Coefficients  of  a Symmetric  Ship,”  J.  Ship  Res., 
Vol.  5,  No.  4,  pp.  1-7  (1962). 
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Since 


by  the  Green  theorem  and 


ff  d * = ff  0i  * 


*ln  = - *Dn  on  So 


from  the  kinematic  boundary  condition,  we  can  show  that 

F‘(C  ' = P ll  [JW°  ^ + (0i  + ^ ^3  6i5  - jjj  ^2  ^ib)  I;]  *1  ds  (?1  > 

So 

Since  we  have  assumed  that  nj  = o(n-,  or  n3),  it  follows  from  Equation  (67)  that 
0|  = Oinj  (=  o(0j,  i=  2,  3,  ....  6).  Hence,  we  let 

F,(e)  = )pu>0JJ n,  0,  ds 

So 

The  above  procedure  for  eliminating  the  diffraction  potential  <pD  from  the  expression  for  the 
wave-exciting  force  and  moment  was  first  shown  by  Haskind36  for  zero  speed  and  was  later 
extended  by  Newman37  for  the  case  of  forward  speeds.  It  is  referred  to  as  the  Haskind- 
Newman  relation. 

Similarly,  we  can  derive 

Fi (m)  = P 21  Sko  Jf  (jwni  + Un35i5-Un2  &i6)0kds 


22  *ko  ff  (</>in  + ^ ^3n  6i5  - ^ *2n  6i6)  K ds 


If  we  express  in  the  form 


Fi(m)=22  ?ko(w2  Aik +JwBik) 


36lfaskind,  M.D..  “The  Exciting  Forces  and  Wetting  of  Ships  in  Waves,”  Isvestia  Akadcmic  Nauk  SSSR, 
Otdelenie  Teknicheskikh  Nauk,  No.  7 (1957):  David  Taylor  Model  Basin  Translation  307  (1062). 

37Newman,  J.N..  “The  Exciting  Forces  on  a Moving  Body  in  Waves,"  J.  Ship  Res..  Vol.  9,  No.  3. 
pp.  190-199  (1965). 
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where  Ajk  are  the  so-called  added  mass  quantities,  and  Bjk  are  the  damping  quantities,  we 
find  from  Equation  (72)  that 

A,k = Rej  ["  5 + § 6)5  ■ ^ 02n  6j6) 0k  r 

o 

B'k  ' lml  ["£  JI  {*'-  + ^ ®3n*i5  - ^ *2n  ««,)  «k  *]  «' 

To  be  consistent  with  earlier  approximation,  we  let 

A i k = B | k = 0 for  k = I,  2 6.  ( ' 

We  can  assume  that  there  is  no  coupling  between  the  motions  on  the  horizontal  and 
vertical  planes  for  slender  ships,  so  that 


A,k  = Bik  = 0 


for  the  following  combinations  of  i and  k 


i = 2,  4,  6 fork  =3,  5 
and 

i = 3,  5 for  k = 2.  4,  6 

As  has  been  seen  previously,  the  hydrodynamic  coefficients  appearing  in  the  equations 

of  motion  can  be  obtained  if  solutions  of  the  velocity  potentials  0j  (i  = 2, 6)  are  known. 

In  the  solution  of  0j.  the  flow  around  each  transverse  section  is  assumed  to  be  two  dimen- 
sional. and  thus  the  variable  x enters  as  a parameter  in  the  expressions  for  0j.  We  then  have 

\k  ■ R“,  [-  A i iy.  *>  * ~ «iS  " A *2u  V.)  K d‘]  i761 

Bi« * lmi [-  Z /dx Ia% ,( A + 8iS - £ 6ie)  o'k  «]  07) 


where  <t>\  is  the  two-dimensional  velocity  potential  function;  the  subscript  N means  the  normal 

derivative  in  the  y-z  plane  only,  / is  the  ship  lengthwise  integral  and  / is  the  contour- 

J L “X(x) 

wise  integral  at  section  x. 
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APPENDIX  B 


APPROXIMATION  OF  DOWNWASH  EFFECTS  ON  TWIN 
SLENDER  BODIES  BY  FORWARD  FINS 


A SWAT  H ship  is  assumed  to  be  replaced  by  the  main  hulls  without  the  vertical  struts 
and  placed  in  a uniform  flow  in  unbounded  fluid  with  a small  angle  of  attack. 

Downwash  effects  by  the  forward  fins  of  a SWATH  ship  will  be  examined  by  assuming 
a horseshoe  vortex  model.  The  vortex  line  eminating  from  the  fin  is  assumed  to  be  located 
at  a distance  of  rr/4  times  the  span  of  the  fin  from  the  body-fin  juncture  at  the  trailing  edge; 
this  is  based  on  the  assumption  of  elliptic  circulation  distribution  along  the  span  of  the  fin. 
This  vortex  line  i-  assumed  to  be  parallel  to  the  body  axis;  the  other  vortex  line  inside  the 
body  is  the  reflected  image  of  the  former  about  a circular  cylinder  which  is  assumed  to  repre- 
sent a main  hull  of  the  SWATH  ship. 

Based  on  these  assumptions,  Pitts  et  a!.28  derived  the  following  two  lift-curve  slopes: 


(<-La’W2  *^W(Bi  ^2 


4?r  (fj  — Tj  ) 


i-  I + ri  ~ rr> 

lC  lo'biv>  <c l-a’w l — -(f]  _ r]  , f] 


where  the  subscripts  T(V)  and  B(V)  respectively  refer  to  the  induced  lift  due  to  the  trailing 
vortices  on  the  tail  fin  and  on  the  body  after  the  forward  fin.  In  these  equations,  i is  the 
tail  interference  factor.  C->  is  the  average  chord  of  the  aft  fin.  f|  is  the  transverse  distance 
from  the  body  axis  to  the  vortax  line  outside  the  body.  T|  and  r-,  are  respectively  the  radii 
of  the  body  at  which  the  forward  and  the  aft  fins  are  attached.  (CIa>  ( and  (CLa i^.,  are  the 
lift-curve  slope  of  the  forward  and  the  aft  fins  as  defined  by  Equation  (42),  and  KW(B)  is  as 
defined  by  Equation  (40). 

For  SWATH  6A.  we  have  (see  Tabic  5)  (CLty)  = 2.18,  (CLa)w^  = 2.0.  KW(B)  = 1.37. 

(A  = 4.48  m.  fj  = 4.24  m,  tj  = 2.29  m.  ^ = 3.01  m.  and  i = — 1 .08  which  is  obtained  from 
Chart  7 in  Pitts  et  al.28  Substitution  of  these  values  into  Equations  (78)  and  (79)  yields 


(C,  qj ) = 0.39 

I.«  B(V) 


Since  the  total  lift  due  to  the  trailing  vortices  from  the  forward  fin  is  obtained28  by 
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L(V)  = |((La'T(V,  +((  i.alB(V)l  3 LI2  A,(f) 


the  additional  lift-curve  slope  to  be  assigned  to  the  aft  fin  is  obtained  by 


1 ( 1 a'(  V l l,(  ‘-a,T(V)  + <C'IalB 


A/f> 


8.05 


= ! , 8 + 0.39)  — = - 0.26  f8. 1 % of  C LQ  2 ) 


which  corresponds  to  less  than  10  percent  of  the  lift-curve  slope  without  the  downwash  effect. 

Another  downwash  effect  which  is  contributed  by  the  Magnus  effect  on  the  body 
induced  by  the  trailing  vortices  and  the  transverse  flow  velocity  due  to  a mutual  blockage 
effect  between  the  two  hulls.  This  effect  will  be  examined  next. 

Ihe  average  circulation  per  unit  span  T on  the  forward  fin  is  obtained  by 

UA/f>cLala 


It  we  denote  the  transverse  flow  velocity  on  one  hull  by  v(x).  then  the  vertical  lift  L and 
pitch  moment  M_  induced  by  the  Magnus  effect  are  given  by 


(81 ) 


(81) 


where  C,  and  xa  < o are  respectively  the  x-coordinates  from  the  center  of  gravity  to  the 
quarter-chord  point  of  the  forward  fin  and  the  rear  end  of  the  body. 

The  transverse  How  velocity  between  twin  bodies  can  be  obtained  from  the  solution  of 
the  problem  for  a slender  body  moving  near  a wall.38  That  is. 


v(  x ) 


UA'lx) 


4ny/b2  - r2(x) 


(83) 


c J.N..  "The  Force  and  Moment  on  a Slender  Body  of  Revolution  Moving  near  a Wall," 
r Model  Basin  Report  2127  (1965). 


74 


Here  b is  one-half  the  distance  between  the  axis  of  the  two  bodies,  r(x)  is 
cross  section  at  x,  and  A(x)  is  the  cross-sectional  area. 

Substitution  of  Equation  (80)  and  (83)  into  (81)  and  (82)  gives 


L = 


r r'(  x ) 

dx 


PU2A,(f)  CLal 
4(f,  -r,) 


where  rj  = r(Cj ),  and 


MP 


x rr1 


dx 


PU2  A,(f,CLal 
4(fl  "V 


xa  b 


pU;A/f)qa, 

4,fl  "V 


- r2(x)  dx 


The  integral 


A 


- - r^(x)  dx  can  be  approximated  for  a slender  body 


V_ 

2b7T 


where  V is  the  volume  of  the  body  after  x = ^ . 
Substitution  of  Equation  (86)  into  (85)  yields 


M 


P 


pU2A,(f)  CLQ, 

4(fi  -v 


the  radius  of  the 


(84) 


(85) 


by 


(86) 


(87) 


A sample  computation  will  now  be  made  for  SWATH  6A.  The  lift  and  moment  generated 
on  the  body  by  the  trailing  vortices  will  be  considered  to  be  borne  by  the  forward  and  aft  fins. 
The  additional  lift  on  the  fins  will  be  represented  as  the  additional  lift-curve  slope  to  each  fin 


pU2A,<f> 

4 (f,  - r,  ) CL“' 


b2  - r 2 ) a 


= — ij2/a  (f)r’’  + A (0  C’  la 

o u ,A1  La  1 LLa2'“ 


p U2  A,(f) 

4 ( fj  - r,  ) ( La  1 


(A,'r> c,  c’l„,  *A2m«2C|„2>  (80. 

Here  subscripts  1 and  2 respectively  indicate  the  forward  and  aft  tins.  C'La  j and  C^a-,  are 
the  additional  lift-curve  slopes  to  be  determined,  and  C2  < o is  the  x-coordinate  of  the  quarter- 
chord  point  of  the  aft  fin.  From  Equations  (38)  and  (89).  j and  C^a-,  are  obtained  by 


Lal  4(fj  - r(  )(C,  - C2) 


b(C,  -«,)  - 1 + V I - 


C'  Al<tlVCLal 

La2  4b7r(f,  - r,  )fC,  -C2)A2(f) 

The  values  necessary  for  the  computation  are  shown  below. 


A,(f)  = 8.05  m2 


A2(f)  = 24.04  m2 


f,  = 4.24  m 


b = 1 1 .43  m 


C,  = 18.27  m 


C2  = 26.75  m 


CLQ,=  4.38/rad 


r,  = 2.29  m 


CLa2  = 3.2/rad 
V = 1020.82  m3 


The  foregoing  values  substituted  into  Equations  (90)  and  (91 ) give 

C'Lal  =-  0.047  (1%  of  Cial) 

CLa2=0n9(3-7%ofCLa2) 

As  can  be  seen  from  the  above,  the  downwash  effects  on  the  body  by  the  trailing 
vortices  are  negligible.  Although  the  results  shown  here  are  applicable  only  to  SWATH  6A 
under  various  assumptions,  it  is  felt  that  the  downwash  effects  by  the  forward  fins  of  a 
SWATH  ship  may  not  be  significant. 
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